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Abstract 

We derive a high energy factorization theorem for inclusive gluon pro- 
duction in A+A collisions. Our factorized formula resums i) all order 
leading logarithms (g 2 h^l/n^))™ of the incoming partons momentum 
fractions, and ii) all contributions {gpi,2) n that are enhanced when the 
color charge densities in the two nuclei are of order of the inverse coupling- 
pi^ ~ g 1 - The resummed inclusive gluon spectrum can be expressed as 
a convolution of gauge invariant distributions W[pi,2] from each of the 
nuclei with the leading order gluon number operator. These distributions 
are shown to satisfy the JIMWLK equation describing the evolution of 
nuclear wavefunctions with rapidity. As a by-product, we demonstrate 
that the JIMWLK Hamiltonian can be derived entirely in terms of re- 
tarded light cone Green's functions without any ambiguities in their pole 
prescriptions. We comment on the implications of our results for under- 
standing the Glasma produced at early times in A+A collisions at collider 
energies. 

Preprint IPI1T-TO8/O68, CERN-PH-TH-2008-074. 



1 Introduction 

Collincar factorization theorems [1] that isolate long distance non-perturbative 
parton distribution functions from perturbatively calculable short distance ma- 
trix elements are central to the predictive power and success of QCD. These the- 
orems can be applied to compute inclusive cross-sections of the form A + B — > 
I(M) + X, where I(M) is a set of heavy particles or jets with invariant mass 
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M and X corresponds to the sum over all possible states (including soft and 
collinear hadrons) that can accompany the object I(M). This cross-section, for 
center of mass energy y/s, can be expressed as [2-6] 



CT AB = / dx * dx b fa/A(Xa,^ 2 )f b / B (xb,H 2 ) 

ab J 



x a XbS \i J \ \M n 

In this equation, f a (b)/A(B){ x a(b)i H >) are the non-perturbative "leading twist" 
parton distribution functions which gives the distribution of a parton a(b) in the 
hadron A(B), as a function of the longitudinal momentum fraction x a ( b ) evolved 
up to the factorization scale p 2 , while the hard scattering matrix element a a b 
can be computed systematically in a perturbative expansion in powers of a s = 
g 2 /47r, where g is the QCD coupling constant. Higher twist contributions to 
this formula are suppressed by powers n of the hard scale M. This factorization 
formula is valid in the Bjorkcn limit when M 2 ~ s >• A 2 CD (where A QCD <~ 
200 MeV is the intrinsic QCD scale). 

Our interest here is instead in a different regime of high energy scatter- 
ing where, for fixed invariant mass M » A QCD , one takes \fs — ► oo and thus 
x a ,b — * 0. We shall call this the Regge-Gribov limit of QCD. An important 
insight is that in this limit the field strengths squared can become very large 
(O (■£-)) corresponding to the saturation of gluon densities [7,8]. The onset of 
saturation is characterized by a saturation scale Q s (x) , which opens a kinematic 
window M 2 ~ Q s 2 3> Aq OD accessible at very high energies. The physics of the 
Regge-Gribov regime is quite different from that of the Bjorken limit discussed 
previously. The typical momenta of partons are ~ Q s 3> A QCD and higher twist 
contributions are not suppressed. These considerations are especially relevant 
for the scattering of large nuclei because the large transverse density of partons 
in the nuclear wavefunctions (proportional to the nuclear radius ~ A 1 / 3 ) pro- 
vides a natural enhancement of the saturation scale, Q 2 (x, A) cx A 1 / 3, . Our goal 
is to derive a formula similar to eq. (1) for inclusive gluon production in the 
Regge-Gribov limit. 

The dynamics of large parton phase space densities in the Regge-Gribov 
limit can be described in the Color Glass Condensate (CGC) effective field 
theory where small x partons in hadrons and nuclei are described by a classical 
field, while the large x partons act as color sources for the classical field [9-11]. 
The lack of dependence of physical observables on the (arbitrary) separation 
between large x color sources and small x dynamical fields is exploited to derive 
a renormalization group (RG) equation, known as the JIMWLK equation [12- 
19]. This equation is a functional RG equation describing the change in the 
statistical distribution of color sources W Y [p] with rapidity Y (= ln(l/a;)). It 
can be expressed as 

d -^ = nw Y [p], (2) 
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where Ti is the JIMWLK Hamiltonian 1 . For a physical observable defined by 
an average over all the source configurations, 

(0) Y =J[Dp]W Y \p]0\p], (3) 

one obtains 

^ = (HO) y . (4) 

We have used here eq. (2) and integrated by parts (using the hermiticity of H). 
The structure of H is such that (HO) Y is an object distinct from (0) Y , so that 
one obtains in principle an infinite hierarchy of evolution equations for operators 
expectation values (0) Y [20]. In the large N c and large A mean- field limit, this 
hierarchy simplifies greatly. When O is the "dipole" operator, corresponding to 
the forward scattering amplitude in deep inelastic scattering, the resulting closed 
evolution equation is known as the Balitsky-Kovchegov (BK) equation [21,22]. 

In refs. [23-25] , we developed a formalism to compute observables related to 
multiparticle production in field theories with strong time dependent sources. 
This formalism is naturally applicable to the CGC description of high energy 
scattering 2 albeit, for simplicity, we considered only a scalar 3 field theory. 
(The corresponding QCD framework was briefly considered in ref. [26].) In 
these papers, the formalism for multiparticle production was developed for a 
fixed distribution of sources, with the assumption that the final results could be 
averaged over, as in cq. (3), with unspecified distributions of sources W Yi [p{\ 
and W Y ^[p2\ (one for each of the projectiles). However, we did not discuss in 
these papers the validity of such a factorization formula. 

In the formalism of refs. [23-25] , one can formally arrange the perturbative 
expansion of an observable like the single inclusive gluon spectrum as 

O [pi, p 2 ] = c ° + Cl -9 2 + C2 -9 4 H » ( 5 ) 

where each term corresponds to a different loop order. Each of the coefficients 
c n is itself an infinite series of terms involving arbitrary orders in (gpi^Y ■ We 
call "Leading Order" the contribution that comes from the first coefficient Co : 

0^,,*] = (6) 

In the case of the single gluon spectrum, the first term co/g 2 has been studied 
extensively. In [24] we developed tools to calculate the next term ci . Following 
this terminology, we denote 

CnloIpi-H = ci , C NNLO [/?i,p 2 ] =c 2 g 2 (7) 



1 The explicit form of this Hamiltonian will be given later in the text. 

2 Although the color sources of each nucleus are independent of the corresponding light-cone 
time, their sum constitutes a time-dependent current. 
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However, this strict loop expansion ignores the fact that large logarithms of the 
momentum fractions X\^, can appear in the higher order coefficients Ci^,... when 
\[s is very large. The term c n can have up to n powers of such logarithms, and 
a more precise representation of these coefficients is 

c n = Vi, In (—) . (8) 

The "Leading Log" terms are defined as those terms that have as many loga- 
rithms as their order in g 2 , 

LLog [ Pl , P2 ] = f>„ 9 2n In" ( J-) . ( 9 ) 

9 n=0 V X 1,2/ 

In this work, we will go significantly further than the Leading Order result, 
and resum the complete series of Leading Log terms. We will prove that, after 
averaging over the sources pi } 2, all the Leading Log corrections are automatically 
resummed by the JIMWLK evolution of the distribution of sources, and that 
the event averaged Leading Log result is given by the factorized expression 

<0> LLog = / [Dpi][Dp 2 ] Wy^-y [Pi] W Yb _ +Y [p 2 ] LO [ Pl ,p 2 ] • (10) 

In this formula, Y is the rapidity at which the gluon is measured, and the 
subscripts Ybcam ± y indicate the amount of rapidity evolution 3 of the source 
distributions of the two projectiles, starting in their respective fragmentation 
regions. 

The expressions M / V bcam ±r[pi,2] in eq. (10) are gauge invariant functionals 
describing the source distributions in each of the nuclei. In analogy to the 
parton distribution functions f a (b)/A(B){Xa(b), P 2 ) we introduced previously, they 
contain non-perturbative information on the distribution of sources at rapidities 
close to the beam rapidities. Just as the latter evolve in /i 2 with the DGLAP [27- 
29] evolution equations, the former, as suggested by eq. 2, obey the JIMWLK 
evolution equation in rapidity which evolves them up to the rapidities Ybcam — Y 
and Ybcam + Y from the nuclei Ai and A 2 respectively. As we will discuss in 
detail, the leading order inclusive gluon spectrum, for given sources p\ t2 , can 
be computed by solving the classical Yang-Mills equations with simple retarded 
boundary conditions. Eq. (10) suggests that the result resumming all the leading 
logarithms of the collision energy can be obtained by averaging over this leading 
order result with the weight functionals W evolved from the beam rapidity to 
the rapidity Y at which the gluon is produced. 

In the Regge-Gribov limit, eq. (10) is the analog of the factorization formula 
eq. (1) proved in the Bjorken limit. While we will prove that eq. (10) holds 
for leading logarithmic contributions at all orders in perturbation theory, we 

3 In terms of the center of mass energy ^/s of the collision (for a nucleon-nucleon pair) and 
the longitudinal momentum components of the measured gluon, one has also - at leading 
log - Ybcam - y = ln( v / s/p + ) and Ybcam + y = ln( v / s/p~). 
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have not attempted to show that it is valid for sub- leading logarithms. There is 
currently an intense activity in computing sub leading logarithmic contributions 
in the high parton density limit [30-36] so an extension of our results beyond 
leading logs is feasible in future. There is another aspect of A+A collisions 
that we have not discussed thus far. Our power counting does not account 
for the so called "secular divergences" [37-39]. These are contributions that 
diverge at least as powers of the time elapsed after the collision. Including 
these contributions will not alter our factorization theorem; it does affect how 
"observables" defined at finite times after the nuclear collisions are related to 
quantities measured in A+A experiments. We will address this issue briefly. A 
fuller treatment requires more work. 

The paper is organized as follows. In section 2, we derive an important for- 
mula for the Next to Leading Order corrections to the inclusive gluon spectrum. 
This formula will play a crucial role later, in disentangling the initial state effects 
from the rest of the collision process. In section 3, we will derive the expressions 
stated in eqs. (2)-(4) for JIMWLK evolution of a single nucleus. Albeit the 
result is well known, our derivation is quite different from those existing in the 
literature [12-19,40-44]. We will obtain our result entirely in terms of retarded 
light-cone Green's functions without any recourse to time-ordered propagators. 
We will show that there are no ambiguities in specifying the pole prescriptions in 
this approach. More importantly, our derivation allows us to straightforwardly 
extend our treatment of the JIMWLK equation to the case of the collision of 
two nuclei. This is discussed separately in section 4 where we show explicitly 
that non-factorizable terms are suppressed and our key result, stated in eq. (10), 
is obtained. In the following section, we will relate our work to previous work 
in this direction and briefly explore some of the connections between the dif- 
ferent approaches. In section 6, we will discuss how one can relate our result 
for the Glasma produced at early times in heavy ion collisions [45,25] and its 
subsequent evolution into the Quark Gluon Plasma. We conclude with a brief 
summary and discussion of open issues. There are three appendices dealing with 
properties of Green's functions in light cone gauge relevant to the discussion in 
the main text of the paper. 

2 NLO corrections to inclusive observables 

Before studying the logarithmic divergences that arise in loop corrections to 
observables, let us derive a formula that expresses the 1-loop corrections to 
inclusive observables in terms of the action of a certain operator acting on the 
same observable at leading order. As we shall see, this formula - albeit quite 
formal - can be used to separate the physics of the initial state from the collision 
itself. 

We have in mind an operator made of elementary color fields, which probes 
multi-gluon correlations. To be specific, for a given source distribution, we shall 
consider the quantum expectation value 

0(x,y) = (A i (x)A^y)) , (11) 



5 



in the limit where the time arguments of the two fields go to +00. We chose 
this particular operator because we wish to study the single gluon spectrum -the 
first moment of the multiplicity distribution- in the collision of two nuclei; it is 
obtained by Fourier transforming this bilinear combination of fields. Note that 
the two fields are not time-ordered. The expectation value of such a product 
can be calculated in the Schwinger-Keldysh formalism [46-48], by considering 
that A l (x) lies on the — branch of the contour and A> (y) on the + branch (A 
representation of the Schwingcr Keldysh contour is shown in fig. 1.) 

c + 

) - 



Figure 1: The closed time path used in the Schwinger-Keldysh formalism. 

This section is organized as follows. We first recall the expression of eq. (f f ) 
at leading order in terms of retarded solutions of the classical Yang Mills equa- 
tions. This result is well known and has been derived in a number of different 
ways. We will then discuss the next-to-leading order computation of this quan- 
tity in the CGC framework. There are two sorts of NLO corrections; these are 
the virtual corrections arising from one-loop corrections to the classical fields 
and the real corrections which are obtained by computing the Q |_ propaga- 
tor of a small fluctuation in light-cone gauge. We will show that C NLO can 
be expressed as a linear operator with real and virtual pieces acting on LO , 
plus an unimportant (as far as the resummation of logs of l/xi i2 is concerned) 
additional term. 

2.1 Leading order result 

We showed in [23] that, at leading order, O is the product of two classical 
solutions of the Yang-Mills equations, with null retarded boundary conditions 4 , 

0^(x,y)=A i (x)A^(y), (12) 

with 

[v v ,t^\ = r, 

lim A^(x) = . (13) 

x ^ — 00 

Here, A denotes the classical field, and J M is the color current corresponding to 
a fixed configuration of the color sources. The current is comprised of one or 

4 The retarded nature of the boundary conditions is intimately related to the inclusiveness 
of the observable under consideration. For instance, if instead of the single inclusive gluon 
spectrum, one wanted to calculate at leading order the probability of producing a fixed number 
of gluons, one would have to solve the classical Yang-Mills equations with boundary conditions 
both at x° = — 00 and at x° = +00 (see [49]). 
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two sources depending on whether we consider only one nucleus or the collision 
of two nuclei - this distinction is not important in this section. It is important 
to note that this current, which has support only on the light-cone, must be 
covariantly conserved, 

[T> M ,J"]=0. (14) 

This means that in general, there is a feed-back of the gauge field on the current 
itself, unless one chooses a gauge condition such that the gauge field does not 
couple to the non-zero components of the current on the light-cone. 

Although one can solve analytically the Yang-Mills equations with these 
boundary conditions in the case of a single nucleus [12,50], this is not possible 
in the case of two nuclei, and one must resort to numerical methods to obtain 
results in this case. Fortunately, as we shall see, the discussion of factorization 
in the case of two nuclei does not require that we know this solution analytically. 

Because the solution of the Yang-Mills equations we need is defined with 
retarded boundary conditions, its value at the points x and y (where the ob- 
servable is measured) is fully determined if we know its value 5 on an initial 
surface E -which is locally space-like 6 - located below the points x and y, as 
illustrated in fig. 2. 

x y 




Figure 2: A locally space- like surface S used to define the initial value of the 
color field. 

Therefore, we will write 

0, o (x,y) = LO [A], (15) 

which means that the observable is considered as a functional of the value of the 
color field on the initial surface E. Note that we use the same symbol for the 
color field and for its initial value on E, although mathematically these objects 
depend on a different number of variables and are therefore different functions. 

5 Since the Yang-Mills equations contain second derivatives with respect to time, one must 
also know the value of the first time derivative of the field on this initial surface. 

6 This means that at every point u£E, the vector n M normal to S at the point u (n^dx/j. = 
for any displacement dx^ on S around the point u) must be time-like. This condition prevents 
a signal emitted at the point u S S, propagating at the speed of light, from encountering again 
the surface X. 
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2.2 Next to leading order corrections 

A detailed discussion of the power counting for moments of the inclusive mul- 
tiplicity distribution can be found in Ref . [23] . The leading order contributions 
to 0(x,y) involves only tree diagrams, which explains why it can be obtained 
from classical solutions of the Yang-Mills equations. As mentioned previously, 
this leading order contribution is of order 0(a~ 1 ) but includes all orders in gp. 
In the rest of this section, we shall study the 1-loop corrections to this quantity, 
that are of order 0(1) in the coupling and to all orders in gp. 

The framework to compute these 1-loop corrections (hereafter called "NLO" ) 
to quantities such as eq. (11) has been developed for a scalar theory in ref. [24]. 
Much of this analysis can be carried over to QCD. To avoid complications such 
as ghost loops, we shall work in a gauge such as the light cone gauge A + = 0. 
Following the discussion for the scalar case, we obtain at NLO, 

NLO (x,y) =A l (x)P(y) + P\x)A>(y)+g l l + {x,y) . (16) 

In this equation, G^ + (x, y) is the — h component of the small fluctuation Schwin- 
ger-Keldysh propagator in the presence of the classical background field A 1 and 
the field f3 l is the one loop correction to the classical field. It is obtained by 
solving the small fluctuation equation of motion 



d 2 U(A) 1 M = 1 d*U(A) 
8A^(x)dA v (x)r" y ' 2dA^{x)dA v {x)dAp{xy ++x 



(17) 

with null retarded boundary conditions : 

lim /3»(x) = . (18) 

x° — > — oo 

Here U{A) is the potential term in the Yang-Mills Lagrangean 7 , obtained by 
writing 

C = £ qU ad - U(A) , (19) 

where £ qua d is defined in eq. (136) of appendix A. We refer the reader to ap- 
pendix A for more details. The source term in this small fluctuation equation 
includes the closed loop formed by the Schwinger-Keldysh propagator Q ++ (x, x) 
to be defined shortly, the third derivative corresponds to the 3-gluon vertex in 
the presence of a background field and 1/2 is a symmetry factor. 

Following [24], we can write the propagator G- + (x, y) in eq. (16) as a bilinear 
combination of small fluctuations of the gauge field whose initial conditions are 
plane waves, 

T^y^r a% Xa (x)a% Xa (y) , (20) 

a, a 



7 Unless one chooses a non-linear gauge condition, U(A) is made of the usual 3-gluon and 
4-gluon couplings. 
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where 



lim a 



±fcA 



d 2 U(A) 
dA^dA^x) 
= e^{k)T a e ±ik ' x 



a±k\a,u{x) = 



(21) 



The sum over A is over the two physical polarizations for the initial plane wave 
and the index a represents the initial color carried by the small fluctuation 
field. In eq. (20), our notation is such that the lower color index (a) represents 
the initial color of the fluctuation, while the upper color index (6 or c) refer 
to its color after it has evolved on top of the classical background field 8 . It is 
important to stress that this decomposition of is valid only if one uses small 
fluctuations that are plane waves in the remote past. Using other solutions of 
the small fluctuation equation of motion (21) would lead to a propagator that 
obeys incorrect boundary conditions. 

The ++ propagator at equal points can be written in a similar fashion as 9 

o«f (*, *) = E / (dpk alb -^ {x) <^ x) ■ (23) 

We note that in a generic gauge, covariant current conservation may require 
the incoming field fluctuation to induce a color precession of the classical current 
J M . This modification of the current will in turn induce an additional contribu- 
tion to the field fluctuation. Our strategy [51-53] to avoid this complication will 
be to perform all intermediate calculations in a gauge where this phenomenon 
does not happen. For instance, on the line x~ — where the color charges mov- 
ing in the +z direction live, on should use a gauge in which A~ = 0. Indeed, 
because the color current only has a + component, covariant conservation is 
trivial in this gauge. A gauge rotation of the final result is then performed to 
return to the light-cone gauge of interest. All effects due to current conservation 
are then taken care of by this final gauge transformation. 



2.3 Rearrangement of the NLO corrections - I 

In this subsection, we will express the small fluctuation propagator Q % l + {x, y) as 
the action of a differential operator on the classical fields A l {x) and A^(y). This 
operator contains functional derivatives with respect to the initial value of the 
color field on S. In the following subsection, we will repeat the exercise for the 
one loop correction to the classical field /3 M (x) and write it in terms of a similar 

8 For future reference, note that quantities with only the lower color index are matrices in 
the adjoint representation of SU(N) defined by 

< fe A>) = <A>)^- (22) 



9 When the two end-points are separated by a time-like interval, there can be an additional 
term contributing to this propagator - see [24] for more general formulas. 
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operator acting on the classical field A^(x). These identities, besides providing 
a transparent derivation of the JIMWLK equation for a single nucleus, will be 
especially powerful in our treatment of nucleus-nucleus collisions. 
Let us begin from the Green's formula for the classical field A 11 , 

= j d 4 y D^ R (x,y) ™ML+BS[A]{x) , (24) 

where Dq P r (x, y) is the free retarded propagator (discussed in appendix A in the 
case of the light-cone gauge) and Bq [A\{x) is the boundary term that contains 
the initial value of the classical field on E. (Boundary terms for the classical and 
small fluctuation fields in light-cone gauge are discussed in detail in appendix 
B.) E + denotes the region of space-time above the surface E. Now, consider an 
operator T (to be defined explicitly later) that acts on the initial value of the 
fields on the surface E, and assume that this operator is linear, which implies 

dAP(y) ~ dAP(y)dA-{y) [V) ' [ ' 

Now apply this operator T to both sides of eq. (24), we get 

T^(x) = J dS D^ R (x,y) ^^S_ TA^y) + TB%[A](x) . (26) 
s+ 

By comparing this equation with the Green's formula for a small fluctuation a M 
(see appendix B), 

«"(*) = / ^y^ x ^ J!^(y) °?<y)+BS[a](x) , (27) 

E+ 

we see that we can identify a^{x) = T A^(x) provided that we have 

B£[a]{x)=TBZ[A]{x). (28) 

Because Bq is a linear functional of the initial value of the color fields on the 
surface E, it is easy to see that the operator T that fulfils this goal is 



J d 3 u [a ■ T u ] , (29) 



where T u is the generator of translations of the initial fields 10 at the point 
it G E. We denote by d 3 u the measure on the surface E (for instance, if E is 
a surface defined by x~ = const, this measure reads d 3 u = du + d 2 u±.) The 
detailed expression of this operator can be obtained by writing explicitly the 



I0p or now, it is sufficient to think of this operator as an operator which is linear in first 
derivatives with respect to the color field on S. 
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Green's formula for the retarded propagation of color fields above the surface 
S, and it usually depends both on the choice of the surface and on the choice of 
the gauge condition. An explicit expression of this operator will be given in the 
next section when the initial surface £ is parallel to the light-cone (u~ = const) 
and when the fields are in the light-cone gauge A + = 0. Therefore, we have 
established the following identity, 



J d 3 u [a ■ T 



0,1* (x) = / d 3 u [a- Tj A»{x) . (30) 



Eq. (30) provides a formal expression of a fluctuation at point x in terms of 
its value on some initial surface £ (in the right hand side of eq. (30), only the 
value of the fluctuation on S appears). This formula is especially useful in 
situations where we can calculate analytically the initial value of the fluctuation 
on E, but were we do not know analytically the classical background field A 
above this surface. 

The single nucleus case is a bit academic in this respect because one can 
analytically compute the background gauge field and the fluctuation at any 
point in space-time. Rather, eq. (30) will prove especially powerful for nuclear 
collisions because in that case one does not have an analytic expression for the 
classical background field after the collision. 

Armed with eq. (30), it is straightforward to write the third term of the right 
hand side of eq. (16) as 

d 3 k 



A, a s 

x [[a_ fcAa • Tj/W] [[a +kXa ■ T v ]A jc (y)] . (31) 



In this equation, the brackets limit the scope of the operators T u „. 



2.4 Rearrangement of the NLO corrections - II 

The terms involving the 1-loop correction (3^ can also be written in terms of the 
operator T u , but this is not as straightforward as for The first step is to 

write down the formal Green's function solution of eq. (17). It is convenient to 
write it as 

(32) 

where [A] (x) is identical to [A] (x) except that all occurrences of the bare 
propagator Dq V r in the latter are replaced in the former by the dressed propaga- 
tor in the background field A^ . This dressed propagator, denoted D^", satisfies 
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the equation 



d 2 U{A) 



D p R Jx,y)^g^S(x~y), (33) 



plus a retarded boundary condition such that it vanishes if x° < y°. 

The second term on the right hand side of eq. (32) is the value /3 would have 
if one turns off the source term (proportional to Q++) in the domain S + above 
the initial surface. It is therefore given by a formula identical to eq. (30), 



J d 3 u [/3 ■ T u ] A" 



(x). 



(34) 



To calculate Pi (x) , let us first make explicit the interactions with the background 
field by writing it as 



d*yD^ R (x,y) 



2 U(A) 



+ 



dA v {y)dAP{y) 
d 3 U(A) 



P p i(y) 



2dA»(y)dAP(y)dA°(y) G ++^ y ^ y \ 



(35) 



This expression is obtained by substituting the expression for the dressed re- 
tarded propagator in terms of the free retarded propagator in the definition of 
/3f. 

Consider now the quantity 



1 f d?k f 

A, a y. 



(36) 



We shall prove that /3f and are identical. Using eq. (30), we can write 

d 3 k 

A, a 



C(«) = ^E J ^ J^mJ r d 3 u[a. kXa .T u ]a% Xa (x) . (37) 



Replace of^ kXa {x) in this equation by the r.h.s of eq. (27). Because the boundary 
term Bq [a + /~\ a ](x) does not depend on the initial value of the classical field A, 
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the action of [a-k\ a • T u ] on this term gives zero. We thus obtain 

d 3 k 

7(2 

A, a ^ 

d 2 U{A) 



A, a s E+ 



y) 



dA-{y)dAP{y) 
d 3 U{A) 



J d 4 yD^ R (x,y) 



dA v {y)dAP{y)dA a {y) 
d 2 U(A) 



[a-kXa ■ T u ]a p +kXa {y) 

[a-kX a -T u ]A' r (y)]a p +hXa {y)} 

C(y) 



8A"(y)dAP{y) 

8 3 U(A) 



+ 



dA lJ (y)dAP(y)dA rT (y) 



(38) 



which is identical to eq. (35). We therefore obtain /3f (x) — (p(x). Combining 
the two contributions /?i and /3 2 , we finally arrive at the compact expression 



d 3 u [0 • T t 



1 r cfik r 
+ 2^J {2*f2E k J [a - feAa ' Tu] [a+fcAa ' 



A»{x) (39) 



We can now use eqs. (31) and (39) to obtain a compact expression for NLO 
corrections to O as 



Cnlo (x, y) 



d 3 u[0-T v 



1 f d 3 k f 
+ 2^J (2v) 3 2E k J ^- hXa ' [ fl+feAa ' T ^ 



+^0 NLO (x,y) , 



(40) 



where we recall that O lo [A] is the same observable at leading order, considered 
as a functional of the value of the gauge fields on the initial surface S. The 
corrective term AO NLO (x, y) is defined by 

A, a 



{ 



a_ fe Aa • T u ]A lh {x)\ [[a +kXa ■ T v ]A jc {y) 

[a +k xa ■ T u ]A tb (x)] [[a_ feAa • T v ]A jc {y)] } . (41) 
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As we shall see later, this term AO NLO does not contain any large logarithm. 
Only the terms in the first and second lines of eq. (40) will be important for our 
later discussion of factorization. 

3 JIMWLK evolution for a single nucleus 

Eq. (40) is central to our study of NLO corrections and of factorization. In the 
rest of this section, we will show how this formula is used to derive the JIMWLK 
evolution equation. In section 4, we will show that it can be generalized to the 
collision of two nuclei. A very convenient choice of initial surface £ in the 
derivation of the JIMWLK equation is the surface defined by x" = e. One 
should choose e so that all the color sources of the nucleus are located in the 
strip < x~ < e. An illustration of the objects involved in eq. (40) and their 
localization in space-time is provided in figure 3. 




Figure 3: NLO corrections in the single nucleus case, seen as an initial value 
problem on the surface x~ = e. The shaded area represents the domain where 
the nuclear color sources live (0 < x~~ < e). The field fluctuations represented 
in red continue to evolve in the region x~ > e until they hit the operator we 
want to evaluate. However, this evolution is entirely hidden in the dependence 
of the classical field upon its initial value at x~ — e, and we do not need to 
consider it explicitly. 



3.1 Gauge choice 

We need first to choose the gauge in which to perform this calculation. Because 
the observable we wish to calculate and everything else in eq. (40) is expressed 
in terms of light cone gauge (A + — 0) quantities, we need to obtain a±k\ a and 
/3 in this gauge as well. However, as previously mentioned, covariant current 
conservation is most easily preserved in a gauge where the field fluctuations 
have no — component. This is because they do not induce a precession of 
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the color current J + while crossing the light cone. We are therefore going to 
adopt the strategy advocated in refs. [51,52,17,19], that consists in performing 
intermediate calculations in a gauge where A~ = and then gauge transforming 
the final result to A + — gauge. 

As discussed in detail in appendix B, if one uses the LC gauge and the 
surface u~ = e as the initial surface, the linear differential operator a ■ T u that 
appears in the identity (30) should be defined as 11 

S t) 
a-T u = 0-(fl(u) o < ( u )) - 7 — — +n(u)o~(u)- 



S(d-{Q(u)A i {u))) y ' x ' 5(£l{u)A~{u)) 

^■W' Wwi)) ' m 

where is the adjoint color matrix 12 that will be defined in eq. (46) . Note 
that this operator in cq. (42) contains a term for each of the field components 
that must be specified on the initial surface to know completely the field above 
this surface. This operator T u can therefore be interpreted as the generator 
of translations of the initial condition for a classical solution of the Yang-Mills 
equations. It is also important to note that the fluctuation field a^{u) that 
multiplies this operator is evaluated just above the initial surface (at u~ = e). 
Therefore, because one does not require its entire history beyond this surface, 
it can in general be calculated analytically. 



3.2 Classical field 

Let us recall the structure of the classical background field itself. As is well 
known, the field in the Lorenz gauge (d^A^ = 0) has no A~ component, and 
therefore fulfills the A~ = condition. Its explicit expression in terms of the 
color source 13 p in given by 

A + {x) = — \-p( x -,x ± ) , A~=A i = Q. (43) 

The gauge transformation that relates the classical background fields in the 
A + = gauge and the corresponding fields in Lorenz gauge is 14 

A" = rtA^n + -rtd^n , (44) 



n We have omitted the color indices in this equation. Q should be understood as a matrix 
in the SU (N) group, and A as a column vector. f2yt is therefore a column vector whose 
components are (QA) C = Q c bAb- 

12 At first sight, CI does not play any role in the definition of T» — the necessity to introduce 
this matrix f2 in the definition of T u is also explained in the appendix B. 

13 Thc density p of color sources is a gauge dependent quantity. When defined in the Lorenz 
gauge, we denote it with a tilde. 

14 In this expression, Q is a matrix in the group SU(N), while A is a matrix in the adjoint 
representation of the algebra SU(N). The product Cl^ACl is a matrix in the SU(N) algebra. 
Note that depending on the context we use the same symbol for an element A of the algebra 
(i.e. a matrix), and for the vector column made of its components Ac on the basis of the 
algebra. The relation between the two is of course A = A C T C . 
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where the tilde denotes fields in the Lorenz gauge; those without a tilde are in 
light cone gauge. Using the light cone gauge condition A + = 0, we get 



d + n = igA + n , 



(45) 



which admits the Wilson line 

Q(x~ , x±) = T exp 



>9 



dz~ A+(z-,x ± )T a 



(46) 



as a solution. Note that because the color sources do not depend on x + , A + 
and depend only on x~ and x±. The solution of the classical equations of 
motion in light cone gauge is then 



A+ =A~ 



, 



A\x) = -Qt(aT ,x±)d i n{x~ ,x±) 



(47) 



We should comment here on the residual gauge freedom of the classical 
solution. The most general solution of eq. (45) is 



ft(x , x^)0(x + , xj_) 



(48) 



where is an arbitrary x -independent gauge transformation. With this more 
general choice, one obtains 

A + = 0, 



a- = -e+a-e 

9 



0t 



(49) 



The arbitrariness in the solution is because the condition A + = does not fix 
completely the gauge and ^"-independent 0's span the residual gauge freedom. 
Requiring that the classical gauge field be of the form given in eq. (47) amounts 
to the choice = 1. This choice is assumed in the rest of this paper. 



3.3 Field fluctuations on the light cone 

To readers familiar with the structure of the JIMWLK Hamiltonian, the struc- 
ture of eq. (40) is already suggestive. In the rest of this section, we will show 
that the leading logarithmic contributions in this formula - terms that are lin- 
ear in the rapidity differences between the projectile and target relative to the 
observed gluon - can be absorbed into a redefinition of the distribution of color 
sources of the nucleus. Our first task towards this conclusion is to compute 
the value of the field fluctuations a±k\ a and (3 just above the light cone on the 
initial surface u~ = e. 
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Let us consider a small fluctuation a M on top of the classical field A^. The 
relation between the two gauges must be modified, 

A" + = + a")fi + ^fit0"fi , (50) 

with 

n = (i + iguj)n, (51) 

where w has components of order unity. Using this ansatz in cq. (50), and 
keeping in mind that A, A ~ 0(g~ 1 ) while a, a ~ 0(1), we obtain the relation 

a^ = f7 t (5^-z ff [w,>]-9^)fi. (52) 

To determine w, as previously, apply the gauge condition a + — 0. This gives 

d + w + ^[w,.4 + ] = 5+ , (53) 
the solution of which can be written as 

Lu(x) = £l(x~ , xj_)f(x + , x±_) + / dz~ Q(x~ , z~; x±) a + (z~, x + , xj_) . (54) 

J — oo 

In this equation / is an arbitrary function that does not depend on x~ , and 
Q.(x~ , z~; x±) is an "incomplete" Wilson line defined by 

0(x~, z~;x±) = T cxp |*5^ dz" ^+(z", a;_L)T a | . (55) 

The arbitrariness in the choice of the function ft, again means that there is a 
residual gauge freedom after we have imposed a + = 0. 

A crucial point in our derivation is how the residual gauge freedom is fixed. 
We need small field fluctuations in order to represent the propagators as in 
eqs. (20) and (23) as bi-linear forms in these fluctuations. These equations are 
valid only if the initial value of the fluctuations a±k\ a are plane waves with on- 
shell momenta; one can check easily that this is true for the free propagators. 
Thus eq. (54) must give plane wave solutions for the field fluctuations in light 
cone gauge when x~ < 0. This is simply achieved by taking plane waves for 
the fluctuation a M in the original gauge and setting the function / to zero 15 . 
Therefore, the requirement that eqs. (20) and (23) be valid leaves no residual 
gauge freedom. 

We only need to know uj on our initial surface S - at x~ = e. Because the 
components of f2 and of a are all of order unity, it is legitimate to neglect the 
values of z~ that are between and e in the integration in eq. (54). For x~ = e 



15 We note that it is also possible to choose a M 's that are not plane waves and a non-zero 
/ to achieve our requirement that a M be a plane wave. This however makes the intermediate 
calculations more tedious. 
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and z < 0, the incomplete Wilson line is equal to the complete Wilson line 
(which has the lower bound at — oo). We therefore obtain 

,o 

u(x = e) = Cl(x±) / dz a + (z ,x + ,x±). (56) 

J — oo 

Note also that when e < x~ , the Wilson line becomes independent of x~ because 
all the color sources are in the strip < x~ < e. This explains why we only 
indicate x± in its list of arguments. 

Once u) has been determined, the — and i components of the fluctuation in 
light cone gauge are determined from those in the A~ = gauge to be 

a 1 = nUti l -d l Lu)n. (57) 



As we shall see shortly when we discuss the leading logarithmic divergences, the 
only quantity we need is 16 

<9 M [fio^fit] = [a" - d^u - ig[u, >]] 

= -0+d-uj - 8 l [h l - d l u\ , (58) 

where we have used eq. (53) and the fact that a~ = in order to eliminate a 
few terms. Using the equation for d + iv, as well as the fact that A + is zero at 
xr = e, we get 

<9 M [Qa^] = 3\uj - d~a+ - d l a l . (59) 

Let us now consider specifically the fluctuations a±k\ a - In the gauge a~ = 0, 
their expression below the light cone reads 17 

~< k Xa^) = ~4(k)T a e ±ik - x , (60) 

with 

r x (k) = 0, 

£ el(fe)et(fe) = - 5 « , 

A=l,2 

et(*) = k -^^ • (61) 



The formulas that govern the light cone crossing in this gauge have been worked 
out in [53]. Using these results, one finds the following expressions for the 
fluctuation fields just above the light cone: 



±ik-x 



l 



«±L» = ^(sJe+W ± (^0 6a (x ± ))— el(fe) e ±ife '* . (62) 



16 Note that (ClACl^) c = Cl c ^Ab, from the definition of the adjoint representation. With the 
notation where A is a column vector, this quantity would also be denoted by (QA) C - 
17 Therefore, a£ 6 fcAo (x) = e^(k)S ab e ±ik - x . 
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Note that for these field fluctuations, one has 



d-a + ±kXa = d^ ±kXa . (63) 

Thus we have 

04 n <fcA« nt ] = d i [d i u-2ai kXa ] . (64) 

Substituting eq. (60) in eq. (56) gives the following expression for u just above 
the light cone, 

uj 

Lu b = T^ba TT~t{(k) e ±ik ' x . (65) 
k ± 

Therefore, 

d^ ±kXa rt] b = -W[e±^af ka e x (k)} , (66) 
where we have introduced the shorthand notation 



/ k k \ 

e[(k) EE ( 5 lrn_ 2 __jzn {k) 



k 

af ka = (6 u -^)n baTl ^-d*n ba . (67) 
3.4 Logarithmic divergences 

Let us recall that our objective is to isolate the leading logarithmic contributions 
to eq. (40). From the structure of this equation, isolating these contributions 
requires that we examine eq. (42) term by term. As we shall see later, the 
contribution in (3 ■ T ( "virtual correction" ) can be derived from the term bilinear 
in T ("real correction"). Therefore, let us concentrate on the bilinear term for 
now. 

To determine the leading logarithmic contributions in the real correction, 
we need to consider the integration over the on-shcll momentum /c M as well. It 
involves an integral 

(68) 



/ 

Jo 



k+ 



which potentially leads to logarithmic singularities both at k + — > and at k + — > 
+oo. Note that wherever k~ appears in the integrand, it should be replaced by 
the on-shell value k~ = k 2 L /2k + . Inspecting the integrand of eq. (40), one sees 
that the k + dependence contains exponential factors 

There is no factor depending on v~ — because the points u and v are both on 
the initial surface S, and thus have equal — co-ordinates. It is clear the integral 
converges at fc+ — ► + thanks to the oscillatory behavior of this exponential. 
On the other hand, when k + — > +oo, the exponential goes to unity and one 
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may have a logarithmic singularity there. However, to truly have a divergence, 
the other factors in the integrand should not have any power of l/k + . 

Let us now examine these. The coefficients in the operator a ■ T u are the 
initial values of £la~, (9 - (f2a*) and 9 M (f2a M ). We need only to keep the coeffi- 
cients that have no power of l/fc+ . One sees readily that this is not the case for 
Qa~ or d~(fla l ) : these two quantities (compare eq. (57) to eqs. (58) and (59)) 
contain a factor k~ <~ l/fc + . 

Thus, as previously anticipated, the only divergence arises when one picks 
up the term d^^a^) both in a ■ T u and a-T v . 

In order to regularize the integral over we must introduce an upper 
bound A + . Physically, this cutoff is related to the division of degrees of freedom 
one operates with in the CGC: the color sources describe the fast partons and 
thus correspond to modes k + > A + , while the fields represent the slow degrees 
of freedom that have a longitudinal momentum k + < A + . Therefore, when one 
performs a calculation in this effective description, the longitudinal momentum 
of all the fields and field fluctuations should not exceed A + , in order not to 
overcount modes that are already represented as part of the color sources p. 
The lower scale in this logarithm is of the order of the longitudinal momentum 
p + of the produced gluon. Therefore, the logarithm resulting from the k + 
integration is a logarithm of A + /p + . 

To pick up the logarithm, one should approximate the exponential by unity. 
This implies that the coefficient of the logarithm is independent of u + and v + or, 
in other words, it is invariant under boosts in the +z direction. As we shall see, 
such perturbations of ^(Sla^ 1 ) can be mapped to a change in the color source p, 
and these logarithms can be absorbed in a redefinition of the distribution W[p\. 



3.5 Real corrections 

Keeping only the term in 9 jU (f2a' i ) in eq. (42), and limiting ourselves to the 
divergent part of the real correction for now, we see that we must evaluate the 
operator 

xE^(^»(^) e * 1 " 11 ) 9, ( aJ *( ,, J e "' !i "' i ) 

a 

r s 2 

x / dn + dv + (70) 

J 5d^(u) bd AZ(u+,u ± ) 5d^(v) ce AZ(v+,v ± ) ■ y ' 

Here, to avoid any confusion, we have written explicitly all the color indices. 
Note also that we have performed the sum over the two polarization states of 
the field fluctuation in this expression 18 . 



18 A useful identity is 



k 
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The object on which this operator acts is the observable calculated at leading 
order, considered as a functional of the initial value of the fields A 1 in light cone 
gauge. In this gauge, the initial values of A + and A~ are zero (provided the 
residual gauge freedom is fixed as explained in section 3.2). Moreover, from the 
set-up of the problem, it turns out that these initial fields do not depend on x + , 

A i (x+,x ± ) = A i (x ± ) , (71) 

and 

^n(«) M ^5(« + , tii) = -d l n( u± ) bd A2u±_) . (72) 

When we restrict ourselves to functional that depend only on x+-independent 
initial fields, we can simply write 19 

du + * = 6 f 73 ) 

sd^(u) bd A%(u+,u ± ) sd*n( u± ) bd A d (u ± ) • 1 ' 



I 



Our goal now is to relate the leading logarithmic contribution we have iden- 
tified to the JIMWLK evolution of the distribution of color sources. As we 
have seen in the previous sections, the initial value of the field in light cone 
gauge has a simple expression when expressed in terms of the sources p or fields 
A + in Lorenz gauge. Therefore, we will try to make the connection with the 
JIMWLK equation in this gauge. To do this, we must relate the functional 
derivative 5 / 5d t n(u±)hdA d (u±) to the functional derivative S/SA + . We begin 
by considering the light cone gauge expression for the classical transverse gauge 
fields given by eqs. (47) and (46). Rewriting A l (x±) more explicitly as 

/■X ~ 
dz~ SJt^.iijtolz-^jJdfz-,^) , (74) 

one observes that a variation 20 5A + (e, x±) of the field in covariant gauge in the 
last x~ bin (of width dx~) leads to a change SA t (x±) of the initial value of the 
gauge field in light cone gauge, given by 

SA i {x±) = -&(x ± )(d l 6A + (e,x ± )dx-^n(x ± ) . (75) 
From this formula, we get the variation of d t n(u±)bdA l d (u±) , 

(5[a^(M ± ) bd ^(M ± )l = -&* ± 6A + {e,x±)dx- . (76) 



19 It is useful to recall that the dimension of a functional derivative operator is Mass _d ^' _D 
where d(A) is the mass dimension of the field with respect to which one is differentiating, and 
D the mass dimension of the space in which this field lives. For instance 

5 o 6 i 
■ Mass"* , —T- ~ Mass 1 . 



8Al(u+,u±) SAl(u±) 

20 It is natural that the size of the bin in which the field A + is changed plays a role here. 
Indeed, because A + is integrated over x~ in the expression of A 1 , a change in a bin of zero 
width produces no change in A 1 ■ Note also that the factor dx~ in cq. (75) is necessary on 
dimensional grounds. 
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Inverting this relation, one obtains 



S f s 

— — — — = - / d 2 x ± G(u± - x±) . (77) 

*a*n(u ± ) M ^(u ± ) J 6A+(e Y ,x ± ) 

Here G(u± — x±) is a two-dimensional propagator whose main properties are 
discussed in appendix C. 

It is important to observe that the functional derivatives on the left and 
right hand side of this equation do not have the same dimensions. This is 
because they are defined with respect to fields that live in spaces with different 
dimensions. On the left hand side, the initial transverse field in light cone gauge 
does not depend on x~ as soon as we are outside the nucleus and is therefore 
a function of u± only. On the right hand side, the field A + depends crucially 
on x~. The e Y argument in the right hand side of eq. (77) is not integrated 
over, and should be chosen as the value of where the last layer of quantum 
evolution has produced its partons. This is the same as the location e of the 
surface E used for the initial conditions, but the subscript Y indicates that it 
may shift as the rapidity Y increases. 

We can now rewrite the operator in cq. (70) as follows 

lln(^) [ dWl/x^(*±,Vx)-~q- -, (78) 

where we have defined 21 

^(xx.Vx)^/ |^ / d 2 u ± d 2 v ± ^a% a (u ± )^l a (v ± ) 

i i 3 3 

ik ± -( U± - v± ) u±-x ± v±-y± (7q) 
(ux-xx) 2 (ux-S/x) 2 ' 1 ' 

From eq. (67), ce^ ka can naturally be broken in two terms. If we keep only the 
first term in each of the a's in eq. (79), we obtain correspondingly 

1 fj2.. J2„. {x\-u\) (y 3 ± -v 3 ± ) 



^)(x±,y ± ) = -g^j ydVLdV 7^ 



u±) 2 (y ± - v_ 



xA'^ui-oJ n(u)W{v) - 1 . (80) 

L J be 

Here the function A y is defined in eq. (167) of appendix C. When we keep the 
first term in the first a and the second term in the second a (or vice versa), we 
get zero because the two terms in a are mutually orthogonal. If we keep the 



1 We performed along the way an integration by parts and used the identity in eq. (164). 
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second term in each of the a's, we obtain 



+ 



d 2 u ± {x\ - «iKj/j - uj) 

(27T) 



1 

8^ 



(x_l - it±) 2 (y ± - it±) 2 
9,{x)Vl\y) - Q,{x)n\u) - n{u)n\y) + 1 

d «ic( t t rj 

xA ij (uj_ -v ± ) \n(u)£tf(v) - 1 . 

L J fcc 



6c 



(81) 



When we add the two contributions, the terms involving A 1 - 7 ' cancel, and we 
are finally left with 



rf 



bc (x ± ,y ± ) = 



d 2 u ± (x\ - u\)(y\ - u\) 
(2-7r) 2 (xj_ - u±_) 2 (y ± - u±) 2 

il(x)ni(y)-n(x)Stf(u)-n(u)rt{y) + l\ . (82) 



This function is precisely the function r] bc (x±,y ± ) that appears in the JIMWLK 
equation [17,19]. 

At this point, a word must be said of the term AO NLO in eq. (40). It is 
given by the difference of two terms that can be obtained from each other by 
exchanging a + k\ a and a^kXa- Going back to the calculation of rj bc (x±, y±), it is 
easy to check that for the calculation of the leading log term these two terms give 
the same result and cancel. Physically this is due charge conjugation symmetry 
- because the classical field is real we obtain the same result by exchanging the 
negative and positive energy asymptotic solutions for the quantum fluctuation, 
AO NLO is the difference between these two and thus cancels out. 



3.6 Virtual corrections 

In the previous subsection, we focused on the real contribution to eq. (40). We 
now turn our attention to the term in j3 ■ T u in eq. (40). Recall that (3^ is 
the one-loop correction to the classical field in the LC gauge and is evaluated 
in eq. (40) at u~~ — e, just above the region occupied by the nuclear sources. 
Mimicking the evaluation of the real contribution, we can write directly 22 



J du + d 2 u± [0 • T u ] = 

U~ —€ 

= J d 2 x ± J d 2 u ± G(x±-u±) 0£(fi(u) M /?£(u) 



6Al(e Y ,x±) 



.(83) 



In (£) v\ Xl _) 



22 Onc can confirm that (P)bd@d an d 9 {^i)hdPd are zero an d therefore cannot appear in 
the operator \J3 ■ T„]. 
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We anticipate that a large logarithm in the k + integral will show up in this quan- 
tity, and we have defined v b (x±) as its coefficient. Note that in this definition 
of the function v b (x±), the value 23 of u~ is u~ = e. 

We begin with the Green's formula for the 1-point function f3^(u), where 
the initial surface is taken at v~ = (instead of v" = e), 

/3"(u) = / d 4 v (u,v) \ „ A 9 y^} , ; P p {v) 
' v 1 J °' fiV ' ; VdA v {v)dAP{v) y y ' 

v~>0 

1 d 3 U(A) 
~2dA v (v)dAP(v)dA' T (v) 



+* »* v ,,» *n"L*„,, S&M ■ ( 84 ) 



By this choice of the initial surface, we do not have a boundary term, because 
f3^ is zero at u~ < 0. The propagator Q^ + {v,v) can be expressed in terms of 
the field fluctuations a±k\ a by using eq. (23). Consider now the Green's formula 
for the fluctuation a + k\ a we introduced in eq. (27), but written this time for an 
initial surface at vT = 0, 



l +k\a( X ) ~ 

y->o 



In this formula, both the fluctuation a + /~\ a and the derivative of the gauge 
potential depend on the background classical field in LC gauge. Let us apply 
to this equation the operator 24 [a-k\ a • T] that substitutes one power of the 
background field by a power of a-k\ a - By defining 

= \ E / -jSm- k J d " v [°-**» • r "] <*a«(«) . (86) 

A < a *V>o 
we obtain for this object the Green's formula 

«.) = / lS,ig.(«.«) [ J^ M yw 



«->0 

1 <9 3 £/(,4) 

' 2 (v) 



oz.(v,v) , (87) 



where we used eq. (23) for the propagator that appears in the source term. We 
see that £ M and /3 M are identical. Therefore, we have proved that 

#(«) = ^ E / / [«-**- • T.] a&J«) . (88) 



23 The value of u + is irrelevant because the 1-point function fi^(u) propagating over an 
x+-independent background field (and with a vanishing initial condition in the past) is inde- 
pendent of u+. 

24 This operator is similar to the operator a ■ T previously denned, but it performs the 
replacement of fields inside the region of the sources, instead of just on the surface of this 
region. 
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Inserting this expression into the definition of v b {x±) 1 we obtain 



A, a _ . f. 

' w >0 

x | d 2 w ± G(x x - ui) Sj(fl(«)M<C(u)) • (89) 

To obtain a divergence at large fc + , we need to tame the oscillations in this 
variable which exist because we have now u~ = e while v~ can be anywhere in 
the range [0, e]. These oscillations are damped only if v~ is in the immediate 
vicinity of u~ = e. As a corollary, note that the left diagram in figure 3 is 
therefore a bit misleading because the tadpole contribution depicted vanishes 
when the upper vertex of the tadpole is below the light cone. In fact, to have a 
leading logarithmic contribution, this vertex of the tadpole must be very close 
to the surface u~ = e, as illustrated in figure 4. 




Figure 4: Leading logarithmic contribution of the tadpole diagram. 
For sufficiently small dx~ , we can use 



lim 

dx~ — ' 



e 

J dv~ [a-kxa ■ %] = a- k xa • T v , (90) 



namely, we recover the operator that substitutes the background field by the 
fluctuation in the last layer at v~ = e. Again, using the eqs. (73) and (77) from 
the previous subsection, we obtain the operator 

J dv + d 2 Vl _ [a_ fc Aa • T„] J d 2 y ± J d 2 v x G(y ± - v±) 

V~ — € 
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When inserted in eq. (89), this gives 

d 3 k 



- ) 3 2E„ 

\,a 



x f d 2 v ± G(y ± -v ± ) dl(n{v) ce a v _ e kXa (v)) 



J d 2 U± _ G(x ± u ± ) d;(n(u) bd a» d kXa (u)) 



5Aj(e Y ,y ± ) 

(92) 



Note that the product of the underlined terms, by themselves are just 



ln^) V bc (x ± ,y ± ). (93) 



The final step in our derivation is to note that when A shares a color index with 
0, we have the identity [17,19,54] 



6At(e Y ,y ± ] 



d:(n(v) ce a»_l Xa (v)) = 0, (94) 



because of the antisymmetry of the adjoint generators of SU (N). We can there- 
fore move the operator S/SAf(€ Y ,y±) immediately after the measure d 2 y ± to 
obtain 

v b (x ± ) = \f d? Vl _ rj bc (x ± , y ± ) , (95) 

J oA c \€ Y , y j_) 

which is identical to the relation between rj bc and v h in the JIMWLK equation. 



3.7 JIMWLK equation 

We shall now combine the real and virtual corrections to write the JIMWLK 
equation. Using the real correction in eq. (78) and the virtual one given by 
eqs. (83) and (95) we can write the total NLO correction, eq. (40), in the form 

°™ = In(^) nO, Q [A + ] (96) 

LLog \P/ 

where we have introduced the JIMWLK Hamiltonian, 

H = \\ d 2 x ± d*y ±71 -^ rj bc (x ± , y ± ) . (97) 

J SAj(e Y ,y ± ) 5Aj{e Y ,x±) 

Although the coupling does not appear explicitly in the Hamiltonian, it is of 
order a s because of the presence of two functional derivatives with respect to 
classical fields that are of order g^ 1 . 
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We noted that the observable O at leading order can be expressed as a 
functional of the classical gauge field A + in covariant jjauge. The average of 
this observable over all the configurations of the field A + , up to NLO, can be 
expressed as 

(O lo + O nlo > = J [DA+] W[A + ] [o LO + O nlo ] . (98) 
At the leading logarithmic level, this can be rewritten as 

(0 LO + NLO ) h =J [DA + ]{[l + AYn]w[A + ]}0 LO [A + ], (99) 

where we denote AY = ln(A + /p+). Note that AF is also the rapidity interval 
between the slowest incoming sources (that have k + <~ A+) and the measured 
gluon. To obtain this equation, one uses the Hermiticity of TL with respect to 
the functional integration over A + . In writing this equation, we have absorbed 
all the leading logarithms of k + into a redefinition of the distribution , 

w[A + ] -» [1 + aywW[.4 + ] . (100) 

This suggests that the distribution W 7 ^" 1 "] should depend on the scale A + that 
separates the modes described as static sources from the modes described as 
dynamical fields in the CGC description. Of course, this is not surprising in an 
effective theory based on such a separation of the degrees of freedom. For this 
reason, it should be denoted as W A+ [A + ]. Therefore eq. (99) can be written as 

<0Lo+0NLo> L = g / [DA + }{[l + \n^^jH]W A+ [A + }}0, o [A + }. (101) 

Because A + is a an unphysical separation scale, the expectation value of ob- 
servables should not depend on this parameter. Differentiating the previous 
equation with respect to A+ and requiring that the r.h.s be zero, we get 25 

^A_iy A+ [I1 = -w A+ [I1. (io2) 

Equivalently, if Y = ln(P + /A + ) denotes the rapidity separation between the 
fragmentation region of the nucleus (located at k + <~ P + ) and the rapidity 
down to which partons are described as static color sources, we have 

±W Y [A + ]=HW Y [A + } 7 (103) 

whichjs the JIMWLK equation that drives the Y dependence of the distribution 
W Y [A+\. 



25 To avoid confusion, recall that H, and hence 8W/dA + , arc of order a B . Therefore, for 
consistency, one should not keep the term proportional to 'H(dW/dA + ) because it is of order 
al and therefore beyond the accuracy of the present calculation. 
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The above considerations also indicate that the distribution W Y [A + ] must 
be evolved to a scale A + comparable to the typical longitudinal momentum in 
the observable to avoid large residual logs contributing to the latter. Therefore, 
at leading logarithmic accuracy, the expectation value of the observable is given 

by 

(0>LL 0g = / [DA + ] W Y [A + ] O lo [.4+] , (104) 

with Y = \n(P + / p + ) the rapidity separation between the beam and the observ- 
able and W Y \A + \ given by the solution of eq. (103). 



3.8 All order resummation of leading logs 

Thus far, we only considered 1-loop corrections that generate one power of the 
large logarithm of P + . On this basis, we deduced an evolution equation for 
W^[.4 + ] using renormalization group arguments. However, the solution of the 
RG equation is equivalent to a resummation of all n-loop diagrams that have n 
powers of large logarithms of p + . We shall here analyze the structure of higher 
loop contributions to confirm whether the all loop resummation performed by 
the RG equation is justified. 

We will not perform here a detailed analysis of these leading n-loop graphs 
to show that we indeed recover the solution of eq. (103). More modestly, we 
will work a posteriori by examining the solution of the JIMWLK equation to 
see what the n-loop graphs that it resums are. Before proceeding, it is useful 
to recall a crucial property of the JIMWLK Hamiltonian defined in eq. (97). 
The operator H contains derivatives with respect to the field A + (e Y ,x±) and 
its coefficients depend on all the fields A + (x~,x±) for < x~ < e Y . For this 
reason, we will denote it Ti(Y), where the endpoint e Y at which the derivatives 
act is related to Y by Y ~ ln(e y ). It is important to note that in a product 
H(yi)H(y2), the derivatives in H(yi) do not act on the coefficients of TL{y2) if 

2/1 > J/2- 

The JIMWLK equation should now be written as 

^W Y [A + ]=H(Y)W Y [A + ], (105) 

and its solution reads 

W Y [A+] =U(Y) W a [A+] , (106) 

with 

(107) 



U(Y) = T Y 



cxp / dyH{y) 
Jo 



In this equation, T Y denotes a "rapidity ordering" such that products of H's 
in the Taylor expansion of the exponential are ordered from left to right in 
order of decreasing y. Wq[A + ] is a non-perturbative initial condition. U(Y) 



28 



is the evolution operator for the Hamiltonian H(Y). Inserting eq. (106) into 
eq. (104), we get 



(O) 



LLog = / [da+]w [a+]uHy) lo [^ 



(108) 



Because Tl(y) is Hermitian, the Hermitian conjugate of the evolution operator 
U{Y) is the same operator with the rapidity ordering reversed : 



U\Y)=T, 



exp / dyH(y) 



(109) 



where T Y denotes the anti-rapidity ordering. The expansion of W to order one 
in TL gives the leading logarithmic one-loop contributions that we have evaluated 
earlier in this section. (See eq. (99) for instance.) 

If one expands it to second order, we see that the leading logarithmic con- 
tributions in the observable at two loops should be given by 



a 



LLog 



dyi 



dy 2 Hfa)Hfa)0, o [A + } 



(110) 



Because y 2 < y%, the derivatives in Ufa) can act on the coefficients rj and v 
of Ufa). Let us first consider the terms where this does not happen, namely 
where the derivatives in Ufa) act directly on C LO [A + ] . These terms correspond 
to the graphs depicted in figure 5. If we look only at what happens below the 
line x~ — e, these contributions are just disconnected products of terms we had 
already at 1-loop. The analysis we performed of the logarithmic contributions 




Figure 5: 2-loop contributions made of products of pieces already encountered at 
1-loop. Although we do not make this distinction in the figure, one of the factors 
is attached at a slightly smaller value of x~ , because the two Hamiltonians in 
eq. (110) are at different rapidities. 

at one loop extends trivially to these terms and it is easy to see that they have 
two powers of the logarithms. 
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In addition, eq. (110) also contains terms in which at least one of the deriva- 
tives in ^(2/2) acts on the coefficients of H(yi). This corresponds to topologies 
of the type displayed in figure 6. Such terms, that have a gluon vertex inside 




Figure 6: Example of term obtained when the derivatives in H(y2) can act 
on the coefficients of Ti(yi). Here, one of the derivatives in H.(y2) acts of the 
function 77 of 7i(yi) and the second derivative in TL{y2) acts on O lo . 

the region where the sources live, have a large logarithm for the same reason 
that the tadpole has a logarithm in the 1-loop terms. Thus one can see that it 
is crucial to properly order the powers of the Hamiltonian TL in rapidity not to 
lose these terms 26 . 

Finally, there also exist at two loops some topologies that never appear in 
cq. (110), such as those of figure 7. The contributions in this figure are 1-loop 




Figure 7: Some of the 2-loop corrections to the observable O that do not appear 
at leading log. 

corrections to the coefficients of the operators T u , v in eq. (40). In other words, 
these terms generate corrections of order a s to the coefficients in the JIMWLK 

2e For instance, if the ordering of the two Hamiltonians in eq. (110) is reversed, we get only 
the terms of figure 5. 
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equation, and do not have double logs of A + . This explains why they arc not 
generated by the leading log formula in eq. (110). 



4 Nucleus-nucleus collisions 

In the previous section, we obtained an expression for resummed leading loga- 
rithmic inclusive gluon observables in a single nucleus in terms of the equivalent 
leading order observable. Along the way, we presented a novel derivation of the 
JIMWLK evolution equation. In this section, we will extend our analysis to the 
case of nuclear collisions. We will show that the leading logarithms of that 
arise in the calculation of loop corrections to the single inclusive gluon spectrum 
can be factored out in the distributions W[.4i~] and Wpl^~] that describe the 
two incoming nuclei. This result will complete a proof of factorization of leading 
logarithms of l/xi,2 f° r inclusive observables in nucleus- nucleus collisions. 



4.1 Inclusive observables at leading order 

As in the single nucleus case, our discussion is valid for an inclusive multi-gluon 
operator O (corresponding to a moment of the multiplicity or energy distribution 
produced in nucleus- nucleus collisions) but for simplicity, we will focus on the 
first moment of the multiplicity distribution - the inclusive gluon spectrum. As 
we discussed in [24,26], the inclusive single particle spectrum in nucleus-nucleus 
collisions can be expressed as 

E P §L = ^^J^ye^H^-i^ + iE,) 

xEW) (A^x)A v (y)) . (Ill) 

A 

Unsurprisingly, the operator (A IJ ,(x)A l/ (y)} is identical to what we considered 
previously in the single nucleus case. In particular, at leading order, the single 
gluon spectrum is evaluated by replacing the two gauge operators in the right 
hand side of the previous equation by classical solutions of the Yang-Mills equa- 
tions. These classical solutions are obtained by imposing retarded boundary 
conditions that vanish in the remote past. The only difference with the previ- 
ous section and with eqs. (13) is that the current J v that drives the solutions of 
the Yang-Mills equations is now comprised of two contributions corresponding 
to each of the nuclei. This is a significant complication in that, unlike the single 
nucleus case, analytical solutions do not exist. However, the classical fields and 
the inclusive spectrum have been computed numerically [55-63]. 

Formally, the single inclusive gluon spectrum at leading order is a functional 
of the LC gauge fields A\ t 2 of the two nuclei on the surface x~ — e and x + = e 
respectively, or of the covariant gauge fields Af 2 in the strips < x~ < e and 
< x + < e (see figure 8), 

E p ^ ^OUAi,A 2 }^0 LO [At,A^}. (112) 
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This quantity does not depend on the rapidity y ~ \n(p + /p ) because of the 
boost invariance of the classical equations of motion [64-66] . 

4.2 One loop corrections 

At 1-loop, eq. (16) can be used again to compute the inclusive spectrum. The 
manipulations in sections 2.3 and 2.4 were not specific to the case of a single 
nucleus. Indeed, we did not specify the detailed content of the current J p in 
section 2. The only requirement for the validity of the final formula is that one 
chooses an initial surface £ which is locally space-like (or light-like at worst). 

We can now exploit this freedom in the choice of S in order to take a surface 
that treats the two nuclei on the same footing. A convenient choice is a surface 
T, with the two branches 

x~ — e , x + < e 

x + =e, x~<e, (113) 
as illustrated by the thick solid line in figure 8. We shall denote the measure 




Figure 8: NLO corrections in the collision of two nuclei. The thick solid line is 
the initial surface where the functions f3^ and a±k\ a are evaluated. The precise 
shape of the small portion of this surface located above the forward light cone 
is not important because its contribution is power suppressed. 

on this initial surface as d£ u . It is simply du + d 2 u± on the first branch and 
du~d 2 u± on the second branch. Similarly, the definition of the operator [a-T u ] 
depends on the branch on which it is evaluated, because the Green's formula 
for the classical fields depends on a different set of initial field components on 
the two branches 27 . It is also important to note that the functional derivatives 
with respect to the initial gauge fields are derivatives with respect to the field 

27 This result is evident from the derivation of the Green's formula in LC gauge discussed 
at length in appendix B. 
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Ai of the first nucleus on the first branch and likewise the field A2 of the second 
nucleus on the second branch. 

We need also to say a few words about the gauge in which the initial fields 
on S are expressed. On the left branch of S (i.e. on the branch u~ = e), we use 
the A + = gauge, while we use the A~ = gauge on the other branch. Using 
different gauge conditions on these two branches is possible because they are 
not causally connected. Similarly, for the propagation of the small fluctuations 
d±k\a and /3, we use the A + = gauge if their endpoint is on the left branch of 
S, and the A~ = gauge if it is on the other side. 

Modulo these obvious changes, eq. (40) is valid in the case of two nuclei and 
we can now express it as 



a 



J dZ u [/3-T u ] 



+ 2 22 J (2ir) 3 2E k J dY ' udY ' v [ a -k^a ■ T u ] [a +kXa ■ T v ] 

A, a E 

+AO NT . n . 



0^ o [Ai,A 2 ] 
(114) 



The first two terms in this formula are illustrated in figure 8. As in the case 
of a single nucleus, the leading logs will cancel in AO NLO because of the charge 
conjugation symmetry discussed previously. 

The leading log piece of the term involving [/3 • TT^ ] can be mapped into the 
corresponding term of the JIMWLK equation in the same way as in the case of 
a single nucleus. Depending on whether we are on the first or second branch of 
the initial surface S, we get two terms which can be expressed together as 



In I J dV u\{ Xl _) — " 



+ ln 



d 2 x ± v%(x±) 



(U^ fe (e y ,x_L) 
S 

SA^ b (e Y ,x±) 



0, o [A+,A^} , (115) 



where v\ 2 { x i-) are respectively the one point functions from the JIMWLK 
Hamiltonian for the two nuclei and likewise, A^ , A 2 are classical gauge fields in 
Lorenz gauge of the first and second nucleus respectively. We have also intro- 
duced a cutoff A~ , that separates the color sources of the second nucleus from 
the dynamical fields. 

There is a subtlety in generalizing the single nucleus derivation to obtain this 
result. In eq. (73), the integration over u + runs from —00 to +00. Now, because 
of the choice of the surface S, this integration runs only from —00 to 0, and we 
must justify that this difference is irrelevant. To simplify the notations in this 
argument, let us use the shorthand a(u + ,u±) = d^{Vl{u)A^{u + ,u±)). In our 
problem, the functional derivative with respect to a(u + ,u±) is only applied to 
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functionate that depend solely on the w + -independcnt mode of a(u + ,u±), 

a(u±) = ^- J du + a(u + ,u±) , (H6) 
where L is the length of the u + interval 28 . When this is the case, we have 

F Hu±)} = y ~T~~TZ. \ F[a(u ± )} . (117) 



Moreover, the result of this differentiation does not depend on the value of u + 
in the l.h.s. Therefore, the subsequent integration over u + merely generates a 
factor L equal to the length of the integration range. We have therefore proven 
that 

/ *+ s^hn 'w-ji - esu F ^ • < 118) 

regardless of the integration range for the variable u + . 

Another possible concern is whether there is a contribution to [(3 ■ T u ] from 
the small portion of the initial surface £ that lies above the forward light cone 
in the region where both u* 1 are positive. It is easy to convince oneself that 
the contribution from this region does not lead to stronger singularities than 
the rest of the initial surface. Furthermore, contributions from this region are 
phase space suppressed due to its small size of order e. 

The leading log contribution of the terms of eq. (114) that are bilinear in 
[a ■ T] is equally simple when the two points u and v belong to the same branch 
of the initial surface S. If this is so, it is straightforward to reproduce what we 
did for a single nucleus, and we find the two separate contributions 



In I ^r) I d 2 x^d 2 y^rii c {x^y^ 



P + JJ ' 5A+ b (e Y ,x ± )SA+ c (e Y ,y ± ) 

0, [A±). 



In ( — ) d 2 x ± d 2 y ± ^ c (x ± ,y ± ) 



P J J " 6A 2b {e Y ,x±)5A 2 (e Y ,y ± ) 



(119) 



Summing cqs. (115) and (119), and expressing v in terms of 77, we obtain the 
leading log 1-loop expression for the single inclusive gluon spectrum to be 



O = 

NLO j , 

LLog 



if ) "I 



0, o [At,A^], (120) 



where Ti\^ are the JIMWLK Hamiltonians of the first and second nucleus re- 
spectively. This equation - assuming we can prove that there are no other terms 
at leading log - is the generalization of eq. (96) to the case of the collision of 
two nuclei. In the next subsection, we will demonstrate that indeed there are 
no other contributions. 



28 Since here this interval is semi-infinite, it is best to consider tt+ e [— L, 0] in all the 
intermediate steps, and to take L — > oo only at the end. 
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4.3 Absence of pre-collision mixings 




Figure 9: Contribution that mixes the two nuclei and may lead to a violation 
of factorization. 

Thus far, we did not discuss the contribution to the bilinear [a ■ T] terms 
where the coordinates u and v belong to different branches of the initial surface. 
This contribution is illustrated in figure 9. If it contains leading log contribu- 
tions, such a term would spoil eq. (120), because it would generate a term that 
mixes derivatives with respect to Af and A^ , thereby precluding any possibility 
of factorization. 

Fortunately, this possibility is not realized because terms where u and v are 
on different branches contain the phases 

e ik+(u~-v~) e ik~[u+-v+) (121) 

in the integral over d 3 k. For generic points u and v in this configuration, neither 
u~ — v~ nor u + — v + are vanishing and these exponentials oscillate rapidly when 
either fc + — > +oo or fc~ — > +oo. Therefore, the integral over k + (or k~) is 
completely finite, and we do not get a large logarithm from this configuration 
of it's and v's. 

The only potential danger might come from the configuration where u or v 
(or both) lie in the small portion of £ above the tip of the light cone. Again, such 
a configuration can at most produce a logarithmic singularity, but is suppressed 
by a small phase space prefactor of order e due to the small size of this region. 
Therefore, eq. (120) contains all the leading log terms that show up in the 1-loop 
corrections to the single inclusive gluon spectrum. 

4.4 Factorization 

Finally, integrating over all the configurations of the nuclear fields Af 2 with 
wcig hts W[Af] and W[A^], and using the fact that the JIMWLK Hamiltonian 
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is Hermitean, we can write the sum of the LO and NLO (leading logs only) of 
the single inclusive gluon spectrum as 

(O lo + O nlo > L = g / [DA+] [DA 2 ] { [l + AY&i] W[A+] } 

x { [l + AY 2 H 2 ] W[A 2 ] } O lo [Af, A 2 ] . (122) 

In this equation, we denote AYi = ln(A^~ /p + ) and AY 2 = ln(A 2 /p~) , where A + 
is the cutoff in the CGC description of the first nucleus, A - of the second nucleus, 
and the longitudinal momentum components of the produced gluon. We can 
now choose the (arbitrary) cutoffs as = p^ and express, as anticipated in 
cq. (10), the leading log part of the NLO result in terms of the LO operator 
convoluted with the appropriately evolved weight functions as 

(0) LLog = J [DAt] [DA2] W Yi [At] W Y2 [A 2 ] 0, o [At,A 2 ] , (123) 

where each of the W[^4 ± ]'s obeys the JIMWLK equation (possibly with different 
initial conditions if the two nuclei are not identical) and Y\ = ln(P 1 + /p + ) and 

y 2 = HP 2 -/ P -). 



5 High energy factorization result in context 

It is useful to consider our result in eq. (123) in the context of related work 
in the high energy limit. Factorization, in the specific sense of our work, was 
proven previously for proton-nucleus collisions in the large N c limit of dipole 
scattering off a large nucleus [67-70]. In the case of nucleus-nucleus collisions, 
there has been recent work by Braun, computing single and double inclusive 
gluon production in a reggeon field theory approach [71]. At present, it is 
unclear how to relate these results to the JIMWLK evolution. A first attempt 
at establishing such a dictionary between cut disconnected diagrams in the 
CGC effective theory and cut Pomerons was discussed in Ref. [23]; see also 
Refs. [72,73]. 

It is important to note that the factorization theorem proven here is valid 
only for inclusive quantities such as moments of the multiplicity or energy dis- 
tributions. In fact, it seems unlikely that these results will extend to discussions 
of total cross-sections and exclusive final states [74-76]. Indeed, it is known 
[24,49] that the retarded nature of the boundary conditions for the fields and 
field fluctuations has a close connection with the inclusiveness of an observable, 
and we have seen in the present paper that the retarded nature of these objects 
plays an essential role in our proof of factorization. Whether the Pomeron loops 
that may play a role in those computations are suppressed for the observables we 
consider is also unclear. Our results certainly suggest that these contributions 
are not important for inclusive moments in nucleus- nucleus collisions, provided 
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the densities pi.2 of color sources are large 29 . 

Another important trend in the literature is computing next-to-leading order 
contributions to high energy evolution. In the reggeon field theory approach [77] , 
we note the very significant work on multi-Regge factorization at NLO by Fadin 
and collaborators [78] which builds on the extension of the BFKL equation to 
NLO [79-81]. In the CGC effective theory, there have been significant recent 
work to include running coupling corrections [30,32-36] culminating in the recent 
NLO extension [31] of the Balitsky-Kovchegov equation. As our result is valid 
for JIMWLK factorization at leading log, these NLO results will be useful in 
attempts to extend our proof of high energy factorization to next-to-leading 
logarithmic accuracy. 

Finally, we should emphasize that JIMWLK factorization proven here is 
far more general and robust in comparison to the k± -factorization often dis- 
cussed in the literature. The latter is concerned with high energy factorization 
at the level of unintegrated k± dependent parton distributions [82-84] and can 
be obtained in the low density limit of JIMWLK factorization [23,85,86]. k±- 
factorization also holds for single inclusive gluon production at leading order 
in proton-nucleus collisions [67,87-91,85]. fc^-factorization was however shown 
to be broken explicitly for quark pair production even at leading order [86] al- 
beit it is restored [23] for large momenta k^^> Q s . Likewise, this breaking of 
factorization is also seen for gluon pair production [68,92]. Though JIMWLK 
factorization remains to be proven for inclusive production of pairs, we antici- 
pate it is far more robust than fcj_-factorization. 

To a large extent, factorization in hadronic collisions is merely a consequence 
of causality : two fast projectiles cannot interact before they collide. Thus the 
objects that describe their content must be universal - independent of the other 
projectile, and of the observable that one is going to measure after the collision. 
However, this general argument does not tell us what information should be 
included in the objects describing the projectiles; indeed, this depends on the 
observable under consideration, and on whether we are in the saturation regime 
or not. In the saturated regime, a given observable will generally be produced 
via the coherent interaction of many partons of the projectiles, which means that 
one will need to know the probability of these multi-parton configuration in the 
wavefunction of the projectiles. In contrast, in the dilute regime, since only 
one parton of each projectile interact, one needs only to know the probabilities 
for 1-parton configurations. This is why JIMWLK factorization is more general 
than /c^-factorization: the distribution W[p] contain enough information 30 to 
calculate the non-integrated gluon distribution, but the converse is certainly 
not true 31 . Similar considerations suggest that JIMWLK factorization may not 
work in the case of exclusive observables. Indeed, inclusive observable usually 

29 If Pi, 2 are not of order g -1 , then the power counting on which our considerations are 
based may be modified. Since it has been argued that Pomeron loops play a role in the dilute 
regime, this leaves open the possibility that these effects may alter our conclusions close to 
the fragmentation region of the projectiles. 

30 It provides information about multiparton correlations such as (p(xi)p(x2) ■ ■ ■ p{x n )\ 
31 Non integrated gluon distributions depend only on 2-parton correlations (p(xi)p{x2)) ■ 
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require less detailed information about the projectiles than exclusive ones 32 . 

The factorization theorem that we have proved here is a necessary first step 
before a full NLO computation of gluon production in the Glasma. Eq. (123) 
includes only the NLO terms that are enhanced by a large logarithm of I/21.2, 
while the complete NLO calculation would also include the non enhanced terms. 
This would be of the same order in a s as the production of quark-antiquark pairs 
[93,94] from the classical field. Note that to be really useful, this complete NLO 
calculation would probably have to be promoted to a Next-to-Leading Log result 
by resumming all the terms in a s (a s ln(l/a;i : 2) n . Now that evolution equations 
in the dense regime are becoming available at NLO, work in this direction is a 
promising prospect. 

6 Factorization, the Glasma and Thermalization 

The Glasma is the non-equilibrium hot and dense matter formed immediately 
in the aftermath of a high energy heavy ion collision [45,95,25]. How this mat- 
ter thcrmalizes is of great importance for a quantitative understanding of the 
phenomenology of heavy ion collisions 33 . We will discuss here the relevance of 
our factorization theorem, present qualitative ideas about its generalization and 
discuss their importance in quantifying the properties of the Glasma. 

At leading order, the Glasma is described by the solution of the Yang-Mills 
equations in the forward light cone with retarded boundary conditions (given 
by the classical fields of the two nuclei before the collision). The produced fields 
have large occupation numbers of order aj 1 and are boost invariant [64,65] . This 
boost invariance of fields implies that the classical dynamics can be described 
by the proper time evolution of gauge fields that live in the transverse plane. An 
interesting consequence of the classical field dynamics is that the chromo-electric 
and magnetic fields are purely longitudinal after the collision [64,45] leading to 
the generation of Chcrn-Simons charge density in the collision [95]. The Glasma 
fields at this order generate only transverse pressure at proper times r > Qs 1 
so it seems impossible that a treatment of the Glasma at this order leads to 
thermalization. 

This is where the small quantum fluctuations of the color field (of order 1, 
compared to the classical field of order g _1 ) become relevant. In an observable 
such as the inclusive gluon spectrum, these quantum fluctuations lead to cor- 
rections that are a s smaller than the leading order classical contribution. As we 
have discussed at length in the previous sections, some contributions of these 
small fluctuations — those that are enhanced by leading powers of \n(\/x\ t 2) — 

32 For instance, in order to study single diffractive processes, one would need "conditional" 
probabilities of multi-parton configurations, where one imposes the condition that no parton 
has been radiated between the rapidity of the projectile to the rapidity where the gap ends. 
This information is not provided by the distributions W[p] that are the basis of JIMWLK 
factorization. 

33 Another important aspect is how jets propagate inside this matter, in order to assess 
issues such as leading parton quenching in jets. 
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can be resummcd and absorbed into universal distributions W[p] that describe 
the high energy evolution of the nuclear wavefunctions. 

But what about the remaining part of these small fluctuation terms, that are 
purely of order a s relative to the classical fields? Our resummation of leading 
logs corresponds to a well controlled approximation provided the coefficients 
d ni in the expansion of eqs. (5) and (8) are truly numbers of order unity. In- 
deed, we have disregarded thus far the terms d n i for i < n, on the basis that 
they do not have as many logs as powers of a s . However, numerical simula- 
tions of the classical Yang-Mills equations with initial conditions that break 
boost invariance show the existence of an instability of the rapidity dependent 
fluctuations [96-98]. In these simulations, it is observed that the small rapid- 
ity dependent perturbations superimposed to the boost invariant classical field 
grow exponentially with the square root of time as 34 



where \i is a quantity of the order of Q s (its precise value depends on the wave- 
length of the fluctuation in the rapidity direction) . This growth has variously 
been interpreted as either a Weibel type [99,97] or Nielsen-Olesen type [100,101] 
instability. The former mechanism in particular has been discussed extensively 
as a possible mechanism for thermalization in heavy ion collisions [102-108]. 
The existence of these unstable modes suggests that our assumption that the 
coefficients d ni for i < n are of order unity is incorrect. 

Our present understanding is that there are three classes among the small 
field fluctuations, that can be organized according to the momentum p v they 
have in the r\ direction : 

• Zero modes (p v = 0) that generate a leading log. That the leading logs 
come solely from zero modes is obvious from the fact that the coefficients 
of the leading logs do not depend on x ± . These terms are already included 
in the resummation we have discussed at length in this paper. 

• Zero modes that do not contribute at leading log because they have an 
extra power of k~ that prevents the divergence when k + — > oo (see the 
discussion in section 3.4). These terms have not been resummed in our 
scheme, and they do not seem to trigger the instability either. They would 
only become relevant in a full NLO calculation, and in resummation of 
Next-to-Leading Log terms [31]. 

• Non zero modes (p v ^ 0). These terms do not contribute large logarithms 
of 1/2:1,2, but they are unstable and grow exponentially as exp( A /7ir). 

It is the latter boost non-invariant terms that are potentially dangerous. While 
also suppressed by a power of a s , they can be enhanced by exponentials of 
the proper time after the collision. Terms that diverge with time are called 

34 The fact that the square root of the proper time, rather than the proper time itself, 
controls the growth of the instability is due to the longitudinal expansion of the system. This 
has also been observed analytically in the study of the Weibel instability [99] . 
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"secular divergences" and some techniques for resumming these divergences are 
well known 35 in the literature [37]. 

Based on the above considerations, let us refine the expansion we wrote in 
eqs (5) and (8), in order to keep track also of powers of exp(^/JFr). We should 
now write ^ 

O [ Pl , P2 ] = -^\c + c ig 2 + c 2 g 4 + ■ ■ ■} , (125) 



with 



n p / 1 \ 

EE^ e(p " )v ^ ln M — • ( 126 ) 

p=i 1=0 v ' ' 

In other words, the coefficients d n i that we have introduced in eq. (8), and 
assumed to be of order unity, are in fact 

n 

dni — ^ ' fnpi , (127) 

p—i 

and can thus grow exponentially in time after the collision. In cq. (126), the sum 
of the number of logs and of factors exp(y/JIf) (this sum is the index p) cannot 
exceed n at n loops. This is because a fluctuation mode cannot be at the same 
time a zero mode (required to generate a log) and a non zero mode (required to 
generate an instability). In this new language, the Leading Log resummation 
that we have performed so far amounts to keep only the term /„„„ in every c n . 

At first sight, one may expect a complete breakdown of the Leading Log 
description when the time 

W~Q. _1 ln 2 Q-) (128) 

is reached. This is the time at which 1-loop corrections become as large as 
the LO contribution. This conclusion can be avoided if one can resum these 
divergent contributions leading to a resummed result that is better behaved for 
r — ► +00. Indeed, it is possible to improve upon the Leading Log approximation, 
by keeping at every loop order all the terms where p = n: this corresponds to 
all the terms where every power of a s is accompanied by either a log or an 
exp(y / 7ir). Thus, let us define 



O 

^LLog+LInst 



00 n / 1 \ 

[p u p 2 ] 9 2n E /«™ e ( "- l)v ^ In* — . (129) 

The subscript "Llnst" is meant for "Leading Instability" . 

In the formalism we have developed in this paper, the growth of small fluc- 
tuations with time can be traced to the action of the linear operator in eq. (114) 
on the classical field. The quantity 

T -**> ~im ~ evW • (130) 



35 Indeed, one can think of the Boltzmann equation as an equation that effectively resums 
a certain class of secular divergences. 
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is a measure of how sensitive the classical field A(x) is to initial condition at 
the point u on the initial surface. If there is an instability, small perturbations 
of the initial conditions lead to exponentially large deviations in the classical 
solutions. We will assume for now that the improved resummation defined in 
cq. (129) can be performed and leads to 

<9 LLog+LI n St = Z[Tu]O^JA] , (131) 

where Z[T U ] is a certain functional of the operator T u . In the r.h.s. we have 
emphasized the dependence of the observable on the initial value of the gauge 
field. This formula can be expressed more intuitively by performing a Laplace 
transform of Z[T„] which reads 

Z[T U ] = J [Da{u)\ e^ d3a Z[a(u)\ . (132) 

Given the structure of a ■ T u in cq. (42), the functional integration [Da(u)\ is 
an integration over the initial fluctuation a M (w) itself and over some of its first 
derivatives. Because T u is the generator of translations of the initial conditions 
on the light cone, the exponential in the previous formula is the translation 
operator itself. When this exponential acts on a functional of the initial classical 
field A, it gives the same functional evaluated with a shifted initial condition 
A + a. Therefore, we can write 

<W LI „ 3t = / [Da(u)] Z[a(u)\ O^JA + a] . (133) 

The effect of the resummation is simply to add fluctuations to the initial condi- 
tions of the classical field, with a distribution that depends on the outcome of 
the resummation 36 . The resummation lifts the limited applicability of the CGC 
approach implied by eq. (128). Indeed, after the resummation, the fluctuation 
a (it) enters only in the initial condition for the full Yang-Mills equations whose 
non-linearities prevent the solution from blowing up. Combining our factoriza- 
tion formula in eq. (123) with the conjectured result of the resummation of the 
leading instabilities, one obtains a generalization of eq. (123) which reads 

(O) LLo g+ Li„ s t - / [DAt] [DA,] W Yi [At] W Y2 [A^] 

x J [Da{u)] Z[a{u)] LO [At +a,A, + a] . (134) 

This formula resums the most singular terms at each order in a s . In comparison 
to the physics of the initial and final state respectively in the collinear factor- 
ization framework, the distributions W[p] are analogous to parton distributions 
while Z[a] plays a role similar to that of a fragmentation function 37 . To prove 

36 In a recent work, using a completely different approach, the spectrum of initial fluctuations 
was found to be Gaussian[109]. 

37 Naturally, this functional has nothing to do with a gluon fragmenting into a hadron. 
Instead, it describes how classical fields become gluons. 
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eq. (134), and to extract the spectrum of fluctuations, one needs to compute the 
behavior of fluctuations on the forward light cone wedge at x T = e, x — > +00. 

Even after the resummations are performed in the initial and final states, 
eq. (134) still suffers from the usual problem of collinear gluon splitting in the 
final state [35]. This however is not a serious concern in heavy ion collisions 
because collinear singularities occur only when one takes the r — > +00 limit. 
In practice, we expect to have switched to a more efficient description like ki- 
netic theory or hydrodynamics long before this becomes a problem. Indeed, the 
initial condition for hydrodynamics, which is specified in terms of the energy- 
momentum tensor T^ v , is an infrared and collinear safe quantity because it 
measures only the density and flow of energy and momentum. It is straightfor- 
ward to re-express our results for multiplicity moments in terms of T^ v . 

A far more challenging problem, that has still not received a satisfactory 
answer, is to understand how the initial particle spectrum - or the local energy 
momentum-tensor - become isotropic and perhaps even thermal. Indeed, a very 
important question is whether this improved resummation, that includes the 
leading unstable terms, hastens the local thermalization of the system formed 
in heavy ion collisions. 

7 Summary and outlook 

In this paper, we have presented a novel derivation of the JIMWLK equation. 
We showed that in this approach the JIMWLK Hamiltonian can be deter- 
mined entirely in terms of retarded propagators with no ambiguities related 
to light cone pole prescriptions. Our approach generalizes easily to the case of 
nucleus- nucleus collisions and we were able to derive the factorization formula 
in eq. (123). This formula is valid to all orders for leading logs in x and to all 
orders in the color charge densities of the nuclei. For this factorization to work, 
it appears crucial to consider an observable that can be expressed in terms of 
retarded fields. Since we had previously linked retarded boundary conditions 
to the inclusiveness of an observable, this emphasizes the importance of inclu- 
siveness for factorization, and the difficulties one may expect when considering 
exclusive observables. 

In view of this, it seems interesting to study whether the factorization theo- 
rem proved here can be extended to less inclusive quantities. One such example 
is the production of two jets that are separated in rapidity by Ay 3> l/a s . 
In particular, can the evolution between the jets be factorized from JIMWLK 
evolution of the wavefunctions as in the case of inclusive gluon production? An- 
swers to these questions will be of great importance in assessing whether the 
early time dynamics in heavy ion collisions leaves an imprint in the long range 
rapidity correlations at later stages. 

We further conjectured the existence of the generalized factorization formula 
in eq. (134). This expression also resums the leading exponentials in time arising 
from the instability of the classical fields to quantum fluctuations on the initial 
light cone surface. The resulting spectrum of fluctuations is very important 
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for determining the subsequent thermalization of the Glasma. Work in this 
direction is in progress. 
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A Gluon propagator in LC gauge 



Consider the QCD Lagrangian to which we add a gauge fixing term proportional 
to (n • A) 2 , 



(135) 



4 -M"- a ■ 2a > 

We are mostly interested in the case where h ■ A = A + , but in fact most of the 
discussion is valid for any vector n M . In order to determine the free propagator 
in this gauge, we need first to isolate the quadratic part of the Lagrangian, 



r ■ — — A a 

i-quad — 2 H 



Ug^ - + -n^n" 
a 



A?, 



(136) 



The free propagator we are looking for is a Green's function of the operator 
in the square brackets. Its calculation is best performed in momentum space, 
where we need to invert 



-g^k 2 + W + -h^n v . 

a. 



(137) 



Because this tensor is symmetric in (/x, u), its inverse must be a linear combina- 
tion of g^ v , k^k v , n^h 1 ' and k^n u +k v n^ . Writing the most general general linear 
combination of these elementary tensors, and multiplying it with eq. (137), we 
finally obtain the following expression for the propagator in momentum space : 



^T(*) = -T3- + 75 



k^k v 



k 2 {n-kf 



a — 



k 2 



+ 



k^i 



k v i 



k 2 (n ■ k) 



(138) 



Note that this expression is still incomplete, because we need to add ze's to the 
denominators in order to make the propagator regular on the real energy axis. 
Doing so amounts to choosing certain boundary conditions for the fields that 
evolve according to this propagator. In this paper, the central object is the 
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retarded propagator, which has all its poles below the real energy axis. This 
amounts to writing: 



k p k v 



k 2 + ik°e (h ■ k + ie) 2 



+ ■ 



k p i 



k u i 



k 2 + ik°e\ ' (k 2 + ik°e){fi-k + ie) 

(139) 

(Our choice for the ie prescription of the n ■ k denominators is indeed retarded 
if n° > 0. We will assume that this is the case.) 

In the case of the light-cone gauge A + = 0, this amounts to choosing a vector 
n p that has h~ = 1 and all its other components zero. Moreover, we work in 
the "strict" light cone gauge, that corresponds to the limit a — > for the gauge 
fixing parameter. The propagator simplifies somewhat in this particular case : 



k 2 +ik e 



k p n u + k v h» 



n • k + ie 



(140) 



Note that this propagator is zero if any of its Lorentz indices is equal to +. 



B Green's formula in LC gauge 

An essential ingredient in our discussion is the Green's formula that expresses 
a field fluctuation in terms of its value on some initial surface. In this ap- 
pendix, this initial surface will be the light-like plane defined by =0, but 
our derivation is more general than that and applies to any initial surface. 

B.l Green's formula for a small fluctuation in the vacuum 

Consider first a small field fluctuation a M propagating in the vacuum. In the 
strict light cone gauge, it obeys 

a+(y) = , 

[V- W]a„( y )=0. (141) 
Recall also that the free propagator D^ R (x, y) obeys 

D% 0tR (x, y) [ by g» v - dffl } = g p »6(x - y) , (142) 

where the arrows indicate that the derivatives act on the left. Now, multiply 
cq. (141) by (x,y) on the left, eq. (142) by a v {y) on the right, integrate y 
over all the domain defined by y~ > 0, and subtract the two equations. One 
obtains 

a p (x)= J A^o, R (^2/)[#;-B y r]«,(!/), (143) 
y ->a 
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where A=A — A- Using the relations 



An B = d»[Ad^ B] , 
A&W B = ^&*[A d V B] +^d u [A d 11 B] , (144) 

we see that the integrand in eq. (143) is a total derivative. Therefore, we can 
rewrite this integral as an integral on the boundary of the integration domain. 
If the derivative we integrate by parts is a d l or d~ , then the corresponding 
boundary is located at infinity in the direction y % or y + respectively. We will 
assume that the field fluctuation under consideration has a compact enough sup- 
port so that these contributions vanish. We are thus left with the terms coming 
from the derivative d + . The contribution from the boundary at y~ = +oo is 
zero, because of our our choice of the retarded prescription for the propagator. 
Therefore, the only contribution is from the boundary at y~ = 0, 

a p (x)= J dy+d 2 y ± D\ , R (x, y) [g^(n- d y ) - \{n>* d v y + n" d £)] a v (y) , 

y-=0 

(145) 

where is a vector such that n ■ A = A (it is the unit vector normal to the 
surface y~ = 0). This formula indicates how the value of the fluctuation at 
the point x is related to its value on an initial surface located at y~ = (Note 
that this dependence is linear since small fluctuations obey a linear equation of 
motion). A priori, it involves the values of all the components of the fluctuation 
on this surface, as well as that of its first derivatives. However, some of this 
information is not necessary because the propagator vanishes when jj, = + and 
because of the gauge condition a + (y) = 0. If one eliminates from the previous 
formula all the terms that are obviously zero and integrate some terms by parts 38 
, we get a p {x) = Bq[o\(x), where Bq[o\(x) is an integral that depends only on 
the value of the field and of some of its derivatives on the initial surface, 

B p [a](x) = J dy + d 2 y ± {[dZD^(x,y)]a-(y) 

y-=o 

-D p ; R (x,y)[d^(y)] -Dg B (x,y) 2d~ a*(y)} . (146) 

Therefore, it appears that in the light-cone gauge A + = 0, and for an initial 
surface x~ = 0, we need to know the initial value of a~ ,d~a l and d^a^ in order 
to fully determine the value of the fluctuation at the point x. This fact is the 
reason why there are only three terms in the definition of the operator T u in 
eq. (42) (but we postpone until the end of this section the explanation of why 
one needs to include the Wilson line in this definition). 



38 The antisymmetric derivatives <9 y and d \ can be eliminated by integration by parts. 

This is not possible for d y since the boundary term does not contain an integral with respect 
to j/~ . This is why we have a term involving the derivative dy D^~ R . 
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Moreover, the first term in the right hand side of eq. (146) can be simplified 
considerably by using the explicit expression of the free propagators in light 
cone gauge : 

Wf>< v) = 5P ~ e ^ - v~) 6 ( x+ - y + ) s ( x ± - v±) • ( 147 ) 

B.2 Green's formula for classical solutions 

There is also a similar Green's formula for retarded classical solutions of the 
Yang-Mills equations. Contrary to the case of small fluctuations, we do not 
assume that the gauge field is small, and we keep all the self-interactions as 
well as the interactions with some external source. Formally, we can write the 
Lagrangian as 

C = Ajuad - U(A) , (148) 

where U(A) is a local polynomial of the gauge field. It contains the 3- and 4- 
gluon couplings and the coupling to the external source. In the A + = gauge, 
the corresponding classical equation of motion is 

[U yg ^-d^ y ]A v { y ) = ^^y (149) 

Then one can follow the same procedure as in the case of small fluctuations, 
and we obtain 

■A p {*)= J d^D^ R (x,y)^^- )+ B p [A]( X ). (150) 

y->0 

Of course, the dependence of the classical field on its initial conditions is no 
longer linear because of the first term in the right hand side; the self interactions 
of the gauge fields lead to an involved dependence on the initial conditions. 

B.3 Green's formula for a M in a background field 

Finally, the Green's formula of eq. (146) can be extended to the situation where 
the fluctuation a^(x) propagates on top of a classical background field *4 M rather 
than the vacuum. The only change is that the free propagator must be replaced 
by the propagator in a background field. The property that its [i = + Lorentz 
component vanishes remains true, because it is a consequence of the choice of 
the gauge. For such a fluctuation, there is also a Green's formula that uses only 
the free gauge propagator, and where the interactions with the background field 
appear explicitly as the additional term 

« P W= / A D^ R {x,y) dA f { y^ {y) a*(y) + £%[«](») ■ (151) 

y->0 
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The derivation of this formula is very similar to that for the classical held . 
We can also rewrite it in a form very similar to eq. (146), i.e. a p (x) = B[a]{x) 
with 

B»{a](x) = J dy+d 2 y ± {[dyD^(x,y)]a-(y) 

-D"-{x,y) [d^(y)] - D*(x,y) 2d~ a l (y)) . (152) 



y-=o 



The boundary term B[a] differs from Bo{a] in the fact that it contains the re- 
tarded propagator D^ v dressed by the background Geld instead of the bare 
retarded propagator Dq^. A crucial difference between the dressed and bare 
propagators is that the simplification of eq. (147) does not occur with the dressed 
propagator. 

In the derivation of the JIMWLK equation, the fluctuations a^{x) one con- 
siders are fluctuations whose initial condition at x° — > — oo are plane waves of 
momentum k. One can calculate explicitly their value on the initial surface, 
which means that we know analytically the quantities a~ , d~a l and d^a^ in 
the r.h.s. of eq. (152). A crucial property is that the initial values of a~ and 
d~a l are suppressed by an extra factor l/fc+, and thus any term containing 
them cannot have a logarithmic divergence when k + — > +oo. This argument 
is correct provided the prefactors of these quantities in eq. (152) do not bring 
factors of k + . There is no problem with the second and third terms, since their 
prefactors is just a propagator. 

However, as we shall see now, the coefficient of the first term can be large 
because it involves the derivative of the propagator. The only case of practical 
interest to us is when the background field above the initial surface is a pure 
gauge field such as the one given in eq. (47). In this particular case, there is a 
simple relationship between the dressed and bare propagators : 

DW{x,y) = n+Oc) D% R {x,y) fi(y) . (153) 

This can be seen by applying a gauge transformation fit to the problem, which 
has the effect of removing the pure gauge background. Using this equation, as 
well as eq. (147), we now obtain 

dj>D»>(x,y) = tf(x) [d»DZ* R (x,v)] n(y) + D»>(x,y) rt(y)d»n( y ) (154) 

The problem is that we take the derivative of the Wilson line ft(y) in a region 
where it is changing very quickly. Only the term with the dy derivative exhibits 
this issue (since the large derivatives are those in the y~ direction), 

D»-(x,y)rt(y)d+n(y). (155) 

From its structure, it is obvious that this term mixes with the second term in 
the r.h.s. of eq. (152) (which, as explained in section 3.4, leads to a logarithmic 
divergence); it would thus be incorrect to keep the latter while not considering 
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the former. There are two ways to deal with this issue: keep track separately of 
these two terms, or try to combine them into a single term. The second option 
is the simplest, and from the above considerations, we know how to achieve it: 
by rotating the fluctuation , — > £!a M , we can rewrite the boundary term as 

BP[a](x) = rt(x) J dy+d 2 y ± {[dy(D^ R (x,y)]n(y)a-(y) 
y-=o 

-D^ R (x,y)[d^(y)a^y)]-D^ R (x,y)2d;Q(y)a i (y)} , (156) 

where we have now only bare propagators. This is why the most convenient 
definition of T u in eq. (42) involves functional derivatives with respect to Oa p 
rather than a M itself 39 . Note that for this discussion to hold, it is only neces- 
sary that the background field is a pure gauge in the vicinity above the initial 
surface, since the derivative is with respect to a coordinate on this initial sur- 
face. Whether the background field is a pure gauge everywhere above the initial 
surface is not important. 



C Two-dimensional free propagator 

In the derivation of the JIMWLK equation, one makes use of several formulas 
involving the bare two-dimensional propagators. These formulas are not new : 
all of them have already been used in one form or another in previous papers 
discussing the JIMWLK equation. We compile them in this appendix, with 
their derivation, as a convenient reference for the reader. 

Let us denote G(x^ — y ± ) a Green's function of the 2-dimensional Laplacian 
operator, 

dlG(x ± -y ± )=S(x ± -y ± ) . (157) 
It admits a simple Fourier representation, 

G(x ± - y± ) = -J*h± e^-(^) _L . (158) 

Note that this object suffers from an infrared problem, which is obvious for 
dimensional reasons: this propagator is a dimensionless object in coordinate 
space, invariant under translations and rotations, and therefore it must be a 
function of fi\x± — y ± \ where [i is some mass scale that was not present in the 
previous equation. 

Derivatives of this propagator do not suffer from this infrared ambiguity. 
Consider for instance 40 

dl G(x± - y± ) = iJ*h± e M.^x) |1 . (159) 

39 Of course, the two ways of defining T„ -with and without the Si— are exactly equivalent. 
But if wc did not include the Si in the definition, the logarithmic divergences would come from 
a combination of the second and third terms of eq. (42), instead of being limited to the third 
term if we include the Si in the definition of T u . 

40 Let us recall that di = -g- = - — . 

x OX; dx z 
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From its symmetries and dimension, it is obvious that this derivative can be 
written as 

diG(x ± -y ± ) = C Xt± ~ V \ 2 , (160) 
{x± - y ± y 

where the prefactor C is dimensionless. Because the derivative of the propagator 
is not infrared singular, the cutoff fj, cannot appear in its expression and C must 
be a pure number (otherwise it would have to be a function of n\x± — y_|_| to 
have the correct dimension). In order to determine the constant, take another 
derivative d l x and integrate over x± the resulting equation over some domain 
£1 of the plane that contains the point y ± . On the left hand side, we get the 
integral of a delta function since G is a Green's function of d\ . We then get 

1 = C [ d*x ± di ( X ^- y \ 2 . (161) 

The right hand side can be transformed by using the 2-dimensional Stokes the- 
orem, leading to an integral on the boundary of f2 (oriented counter-clockwise) 

l-C / (162) 
Jan (x±-y ± y 

where e y is completely antisymmetric (e 12 = 1). The contour integral in this 
equation is a topological quantity, that depends only on the winding number of 
the contour 9f2 around the point y ± . Thus, it is best calculated by deforming 
dfl into the unit circle around the point y ± . We get easily 

1 = 2nC . (163) 

Thus we have 

<^-»i)-s^^?- (164) 

The second derivative of the propagator is also useful in the derivation of 
the JIMWLK equation. By applying d 3 x to the previous equation, one obtains 



didlG(x ± -y ± ) = 1-di '' " 



2tt x {x^-y^Y 
1 



2n(x±_ - y ± ) 2 



gii _ 2 ( x ± -y l ±)( x i -v±) 

(x± - y ± ) 2 



(165) 



This formula, although perfectly correct for x± ^ y ± , is incorrect at the point 
x± =y±- In order to see this, take the trace over the indices i and j. In the left 
hand side, we have the Laplacian of the propagator, i.e. 5(x± — y ± ), while the 
right hand side would give zero. Thus the full formula for the second derivative 
is 

dj c diG(x ± -y ± ) = S -^S(x ± -y ± ) + ^-A ij (x± - y ± ) , (166) 
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with 



W{x ± -y ± ) = 



1 



(x±~y ± ) 2 



S ij -2 



(g±-»i)(gj-yi) 



(167) 



This function A 4 - 7 obeys an interesting identity. By integration by parts, one 
can check that 

d 2 u ± d 2 Vl _ (x i L -u i L )(y 3 L -v 3 ± ) tfrf r( . 
1 f fu± {x\-u^)(y\-u\) 

(271)2 7 ( 27r )2 ( :E± _ M± )2(y ± _ M± )2 



(27r)2(27T) 2 (x ± -u ± ) 2 (y ± -^) 2 
Using now eq. (166), we obtain the following identity, 



S ij 6(u± - v ± ) .(168) 



/ 



d 2 u± d 2 v± (x\ - «^)(y^ - v j j_) 



(27r) 2 (27r) 2 (x ± -u ± ) 2 (y ± -v ± y 



8 %3 1 
— <5 (u ± - v ± ) - — A y (tt ± - w ± ; 

Z Z7T 



= . 
(169) 

Let us also provide an alternate representation of the 2-dimensional propa- 
gator that is sometimes helpful. Let us start with the integral 



rf 2 «i u\ - x\ u\ - y\ 
(2tt) 2 (uj_ - x±) 2 {u ± - y ± f 



-I 



d u± 



d l u G(u ± -x ± ) d:G(u ± - y± 



(170) 

The integral in the right hand side can be performed by parts, since it leads to 
the Laplacian of a propagator, which is a delta function. Thus, we obtain the 
identity 

G{x± - y_J - - / — — - — . (171) 

J (2tt) 2 (u± - x ± ) 2 (uj_ - y ± ) 2 

Note that the integral over u± suffers from the same infrared problems that we 
have already mentioned at the beginning of this appendix. 
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